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ABSTRACT 

The main aims of this research are to extend the averaged Lagrangian 
method Of describing small signal wave propagation and nonlinear wave 
interaction, developed by earlier workers for cold plasmas, to the more 
general conditions of warm collisionless plasmas, and to demonstrate 
particularly the effectiveness of the method in analyzing wave-wave 
interactions. The tneory is developed for both the microscopic descrip- 
tion and the hydrodynamic approximation to plasma behavior. First, a 
microscopic Lagrangian is formulated rigorously, and expanded in terms 
of perturbations about equilibrium. Two methods are then described for 
deriving a hydrodynamic Lagrangian. In the first of these, the Lagrangian 
is obtained by velocity integration of the exact microscopic Lagrangian. 

In the second, the expanded hydrodynamic Lagrangian is obtained directly 
from the expanded microscopic Lagrangian. As applications of the micro- 
scopic Lagrangian, the small-signal dispersion relations and the coupled 
mode equations are derived for all possible waves in a warm, infinite, 
weakly inhomogeneous magnetoplasma, and their interactions are examined 
for propagation nearly parallel to the static magnetic field, and exactly 
perpendicular. As examples of the use of the hydrodynamic Lagrangian, 
the coupled mode equations are derived for interactions among two elec- 
tron plasma waves and one ion acoustic wave; among one electron plasma 
wave, one ion acoustic wave and one ordinary wave, and among two electron 
plasma waves and one Alfven wave. The work concludes with a brief 
discussion of specific cases of wave-wave interaction most likely to 
be worth subjecting to detailed numerical analysis, and some suggestions 
for extension of the Lagrangian method to bounded, inhomogeneous, weakly 
turbulent plasmas of practical interest. 
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Chapter I 


. INTRODUCTION 

1.1 Motivation of the Research 

During the last forty years, the theory of plasma waves has 
developed from the simple analysis of Tonks and Langmuir of electron 
plasma oscillations to very comprehensive treatments based on kinetic 
theory. A wide variety of modes have been predicted by linear theory 
for hot and cold magnetoplasmas, and have been studied experimentally, 
particularly during the last decade. The stage has now been reached 
where it may reasonably be stated that, at least for stable modes in 

homogeneous plasmas, the relevant dispersion relations have been very 

2 

well verified. 

The linear approximation is, of course, inadequate to describe 
unstable waves: the predicted exponential growth of small perturbations 

must ultimately be limited by nonlinear effects. As a result, nonlinear 
and turbulent processes are found to play an important role in many 
laboratory and space plasma situations of considerable interest. For 
example, turbulent heating and anomalous diffusion are vital processes 
in the heating and confinement of plasmas, and the possible achievement 
of controlled thermonuclear fusion conditions. These, together with 
such phenomena as charged particle diffusion out of the van Allen belt, 
and the dynamics of the earth’s bow shock, are highly dependent on 
nonlinear effects and illustrate the importance of understanding them. 

A further example is the need to understand propagation characteristics 
of electromagnetic waves through turbulent plasmas, either from the point 


1 



of view of communications or diagnostics. The description of nonlinear 
effects is therefore one of the central problems in modern plasma physics. 

Unf ortunately, nonlinear phenomena lead very quickly to formidable 

3-5 

theoretical and algebraic complexities. Refuge may be sought in 
computer simulation of the problems, but this can be extremely expensive 
if complicated situations are to be modeled realistically, and it is 
easy for the physics to be obscured. There is consequently very strong 
motivation to elucidate the mechanisms involved with improved analytical 
tools. What is required is a formalism which will simplify the analysis 
to the greatest possible degree. Such a formalism is provided by the 
Lagrangian methods which form the subject of the present work. 

1 .2 Review of Nonlinear Wave-Wave Interaction 

If two waves are propagating as exp[ i(°^t-Jk ’r) ] and 
exp[i(o3 t-k^-^r)], nonlinearity of the plasma will manifest itself in 
the generation of beat waves propagating as exp[i(d3 t-k *r)], where 

3 3 

0 ) =03 ± 0 ) and k =k ±k . If a beat wave happens to be a normal mode of 
5 i 2 ~~j — l ""fc! 

the plasma, synchronism is said to occur, and an exchange of energy and 
momentum will take place among the three waves. The present work is 
concentrated on this nonlinear wave-wave interaction process. It is 
not the only significant nonlinear process; the beat wave can also 
interact with charged particles moving at its phase velocity 03 /k . 

3 ~3 

£ 

Nonlinear wave-particle interaction then occurs. 

A simple degenerate form of nonlinear wave-wave interaction was 

7 

explored by Pierce in connection with the traveling wave tube. In 

g 

1956, Manley and Rowe established an important relationship between 
the powers at different frequencies flowing into a nonlinear capacitor 
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which is now known as the Manley-Rowe relation. Its generalization 
to waves is extremely useful in the analysis of wave-wave interaction 
because it implies certain relations between the coupling coefficients 
in the equations describing the wave coupling. ^ By 1960, amplification 
due to wave-wave interaction, so-called parametric amplification, had 

10 

already been studied extensively in the context of microwave devices. 

Scarcely a paper on parametric amplifiers appeared without mentioning 

the relevant Manley-Rowe relations. Since then, parametric amplification 

has been studied in other fields which involve wave interactions in 

11 12 

nonlinear media, e.g., nonlinear optics and microwave acoustics. 

Nonlinear wave-wave interaction in plasmas seems to have been 

13 

considered first by Sturrock , in 1957- He obtained Manley-Rowe type 

energy transfer relations among four Langmuir waves, and calculated 

nonlinear corrections to the frequencies of the waves. Subsequently, 
lip 

Kino analyzed the nonlinear coupling between three waves propagating 

on a cold plasma column, with and without an infinite static axial 

15 

magnetic field. Sturrock later showed that the Manley-Rowe relations 

represent a special case of the action-transfer relations applicable 

to any system describable by a Hamiltonian. Petschek pointed out the 

possibility of having instabilities due to nonlinear interactions in 

systems stable to small perturbations. The nonlinear interaction of 

longitudinal electron and ion oscillations was studied by Oraevskii and 
17 

Sagdeev . Following these pioneer works, a large number of papers have 
appeared. Many of them will be mentioned in succeeding sections. 

Depending on the situation of interest, two theoretical approxi- 
mations have been made when analyzing wave-wave interaction: the coherent 


3 



phase approximation and the random phase approximation. The former 
is based on the assumption that each of the interacting waves is 
coherent . In plasmas, many waves are available to satisfy the synchronism 
conditions, and the phases of the waves of interest lose coherence due 
to wave-wave interaction with other existing waves, wave-particle inter- 
action or other higher order processes. However, if significant energy 
transfer among the three waves can occur in a time small compared with 
the time required for changes in phases, the waves may be considered as 
coherent, and the coherent phase approximation is adequate. On the other 
hand, if the correlation among the phases of the waves vanishes due to 
the random phase change before appreciable transfer of energy from one 
wave to another, then the modulus of the wave amplitudes averaged over 
the phases is of interest, and the random phase approximation should be 

18 

used. This situation was first considered by Galeev and Karpman, and 

5 19 

has been used extensively in turbulence theory. ' The present work 
is treated by the coherent phase approximation. 

The most important equations to be obtained in the analysis of 
wave-wave interaction are the coupled mode equations, from which spatial 
and temporal variations of amplitudes of the interacting waves can be 
calculated. The Manley-Rowe relations follow from them, and are useful 
in understanding certain features of the parametric interaction. It 
also provide a check on the tedious algebraic manipulations usually 
required in the derivation of the coupled mode equations. We shall now 
discuss several of the mathematical methods that have been proposed and 
used to obtain them. 


4 



1.2.1 The Iterative Method 


The most direct line of attack on wave-wave interactions is 

the iterative method. The procedure is to expand up to second order 

\ 

in the perturbation amplitudes the current equation and the equations 
describing the particle dynamics, i.e., the single particle equation of 
motion for cold plasmas, the moment equations for the hydrodynamic 
approximation, or the Vlasov equation for the microscopic treatment of 
warm plasmas, and to simplify second order terms by the use of small 
signal relations and the synchronism conditions. Combination of the 
resulting equations with Maxwell's equations gives the coupled mode 


equations and the Manley-Rowe relations. 

Up to now, the majority of analyses of wave-wave interaction have 

been done by this method, partly because it is sure and well known 

from the fields of microwaves and optics, and partly because other 

methods have not been well developed. Specific cases treated by either 

the cold plasma or the hydrodynamic approximation are: interaction among 

17,20 

three longitudinal waves, ’ among two ordinary waves and one Langmuir 

21 22 

wave, among two ordinary waves and one extraordinary wave, among 

23 24 25 

three modified ordinary waves, ’ among three extraordinary waves, 

26 

among two circularly polarized waves and one Langmuir wave, and among 

27 

three right-hand polarized waves. The theory has been extended to 

the microscopic treatment for the specific cases of: interaction among 

28 

three longitudinal waves, among two transverse waves and one ion 

29 30 

acoustic wave, y among two ordinary waves and an extraordinary wave/ 

31 32 33 

and among three cyclotron harmonic waves. ' Stenflo has considered 
generally the interactions of waves propagating parallel to the static 


5 




magnetic field, and the extension to waves with their wave vectors 

slightly oblique to the magnetic field have been investigated by Kim, 

3k 

Harker and Crawford. 

The iterative method has the advantage of being straightforward, 
and can be applied to systems for which the Lagrangian approach to be 
described below does not apply, i.e., a Lagrangian function is 
not derivable. In practice, however, it is extremely tedious. For 
warm plasmas, especially, the charged particle velocity distribution 
function has to be expanded to second order in perturbation by solving 
the differential equations derived from the Vlasov equation, with a 
driving term containing the first order distribution function. Since 
the driving term is already complicated, the entire process requires 
very laborious algebraic manipulation. 

1.2.2 The Averaged Lagrangian Method 

In this subsection, we shall discuss the averaged Lagrangian 
method and its application to cold plasmas. The warm plasma case has 
formed the main topic of our work, and will be discussed separately in 

35 

Section 1 . 3 . The procedure is as follows. First, the plasma pertur- 
bation parameters, i.e., the position vectors of particles and the 
fields, are expanded in terms of a sum of sinusoidal perturbations from 
equilibrium whose amplitudes, frequencies, and wave vectors are assumed 
to vary slowly in space and time due to the nonlinearity. The Lagrangian 
is then expanded in terms of these perturbations, and averaged over 
space and time so as to remove rapidly varying terms. The Euler-Lagrange 
equations derived from the zeroth order Lagrangian yield the dynamical 
equations for the equilibrium state. The first order Lagrangian vanishes. 
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Variations of the second order Lagrangian with respect to the amplitudes 
and phases give the small signal equations and the equation of action 
conservation in each mode. Similar variations of the third order 
Lagrangian give the coupled mode equations and the Manley-Rowe relations. 

For the method to be applicable, it must be possible to set up a 
Lagrangian density for a given system, i.e., such that application of 
Hamilton' s principle leads to a set of Euler-Lagrange equations corres- 
ponding to those chosen to describe the plasma behavior. It is well 
known that the description of charged particles in an electromagnetic 
field satisfies this requirement. Sturrock^ used Hamilton's principle 
for electron beams and for cold plasmas to show that the first order 
Lagrangian vanished to within a total divergence, and that the second 
order Lagrangian gave the required small signal equations. 

The averaging process has been discussed by several authors. 

70 

SturrocK introduced a parameter in which the Lagrangian and the dynamical 
variables are cyclic, and showed that the time averaged Euler-Lagrange 
equation obtained from variation of this parameter yielded the action 

39 

transfer equations. Galloway considered the equation for energy con- 
servation averaged in time and space in order to get the wave coupling 

40 

coefficient. Whitham studied a Lagrangian averaged over the fast 
varying waves in a slowly varying medium. The dynamical equations were 
obtained by variation of the averaged Lagrangian with respect to the 

wave amplitudes. From the variation in the phase of waves he obtained 

4-1 

an adiabatically conserved quantity. Bretherton and Garrett showed 
that this adiabatically conserved quantity was the wave energy density 
divided by frequency. This conservation equation is precisely the 
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38 

action-transfer relation previously obtained by Sturrock , who had 
concluded that the Manley-Rowe relations are a special case of the action- 

42 

transfer relation. Dewar has since shown rigorously that this conser- 
vation relation holds to all orders in e , where the small parameter 
e is the ratio of the period (or wavelength) of the fast varying wave 
to that of the slow variation. 

43 

The averaged Lagrangian method has been applied by Dougherty to 
obtain equations for ray tracing, to demonstrate conservation of wave 
action, and to derive nonlinear coupled mode equations in a relativisti- 
cally covariant form. While Dougherty presented the application of the 
method in general terms, specific cases have been studied by various 

44 

authors. Galloway and Crawford have applied it to the nonlinear inter- 
action of waves at arbitrary propagation angles in an infinite magneto- 
plasma, and illustrated its effectiveness in obtaining explicit 

35 

expressions for the coupling coefficients. Dysthe not only gave a 
clear presentation on how the averaged Lagrangian method is applied to 
nonlinear wave interactions, but also studied the self-action effects 
and the decay problem of a finite amplitude wave into sidebands by this 

method. The coupled mode equations for Langmuir and ion acoustic waves 

45 

have been derived by Vedenov and Rudakov. Dougherty considered the 

46 

waves in a cold unmagnetized plasma as an example of his earlier study. 

It is clear from the work cited so far that the Lagrangian method 
is very efficient for treating nonlinear wave-wave interaction. Its 
- advantages over the iterative method derive mainly from the following 
features. First, in the iterative method, the coupled mode equations 
are obtained as the end-product of tedious algebraic manipulation. 


43 
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Furthermore, the Manley-Rowe relations are obtained only from the coupled 
mode equations. In the averaged Lagrangian method, however, general 
expressions for the coupled mode equations and Manley-Rowe relations 
appear naturally before specifying an explicit expression for the 
Lagrangian function for a particular system. Further, the third order 
Lagrangian itself gives the coupling coefficient, and the expansion of 
the Lagrangian only needs relatively simple and straightforward algebra. 
1.2.3 The Hamiltonian Method 

As the Manley-Rowe relations indicate, wave energy appears 
as an important quantity in nonlinear wave interaction. Therefore, the 
development of a Hamiltonian method is quite natural . This approach was 
initiated by Sturrock 1 ^' ^ in his study of the nature of the Manley-Rowe 
relation attributable to any system which may be described by a Hamiltonian. 
The application to wave-wave interaction is similar to that of the 
Lagrangian method. The Hamiltonian generates higher order perturbation 
terms which are averaged over the fast varying temporal and spatial 
variations. The canonical equations from the second order Hamiltonian 
give the small signal equations, and the equations describing the non- 
linear effect are obtained from the third or higher order Hamiltonian. 

In this method, there are two approaches dependent on the choice 

of canonical variables: the first uses action-angle variables, and 

is 17 

was initiated by Sturrock. ’ The Euler-Lagrange equation obtained 
by variation of the Lagrangian with respect to angle is the action 
transfer relation. The other canonical equation gives the frequency- 

shift relation. Even though this frequency-shift relation is obtainable 

35 

by the Lagrangian method, it emerges more naturally in the Hamiltonian 
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method. This approach is analogous to the averaged Lagrangian method, 
but does not provide the coupled mode equation obtainable by the 

49 

Lagrangian method. In order to overcome this disadvantage, Harker has 
recently developed another approach more closely parallel to the second 
quantization method (see Section 1 . 2 . 4 .) than the Lagrangian method. 

This second approach uses for canonical variables the classical para- 
meters corresponding to the probability amplitudes in quantum mechanics. 
These parameters and their complex conjugates are then canonical conju- 
gates, and the canonical equations yield the coupled mode equations. 

He has applied this approach to several illustrative cases: the inter- 

actions of two longitudinal waves and one transverse wave, of three 
longitudinal waves on a positive column, and of three longitudinal waves. 
1 . 2.4 The Second Quantization Approach 

Although the concepts of quantum mechanics, in particular 

18 45 50-52 

the idea of quasiparticles, have been used by Soviet scientists, ’ ’ 

it was not until recently that a formal quantum mechanical method was 

53 

developed for wave-wave or wave-particle interaction problems. The 
procedure consists of expanding the Hamiltonian about equilibrium, and 

54 

then quantizing the fields in order to express the Hamiltonian in 
quantum mechanical form. The transition probability from one state to 
another can then be obtained by the Fermi golden rule in the time 
dependent perturbation technique. For a classical plasma, the transition 
probability is taken to the classical limit. 

55 

The method was originally developed by Pines and Schrieffer. 

Walters and Harri s 56 used the method to study the nonlinear interactions 
of longitudinal waves in a magnetized plasma. The static magnetic field 
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was assumed to be so strong that the electron motion was essentially 
one-dimensional. It was demonstrated that the classical limit of three- 
wave interaction among plasma waves yielded an earlier result of Aamodt 
and Drummond . ^ Krishan and Selim studied the interactions of trans- 
verse waves^ and of transverse and longitudinal waves'^ in an 

unmagnetized cold plasma. Three-and four-wave interactions in a cold 

60 

plasma with no magnetic field have been analyzed by Zakharov. 

The advantages of the method are that, if the Hamiltonian is obtained 
in the quantized quantities, obtaining the transition probability is in 
principle straightforward. While the classical treatment is usually 
confronted with the conflicting choice of coordinate systems, i.e., 
Eulerian coordinates for fields and Lagrangian coordinates for particles, 
such complications from the coordinate systems do not occur in this 
method. The classical limiting procedure is sometimes extremely tedious, 
however. Furthermore, unnecessary quantum mechanical steps have to 
be followed even for problems which can easily be solved by classical 
methods. If there are only a few monochromatic waves in a plasma, then 
the use of the Fermi golden rule to calculate transition probabilities 
is no longer justified. Therefore, this method is primarily for the 
case where the random phase approximation has to be made. 

1 .3 Review of Lagrangian Method for Warm Plasmas 

For nonlinear wave interactions in cold plasmas, the efficiency of 
the averaged Lagrangian method has been demonstrated by many authors. 

The algebraic manipulation in the iterative method is much more tedious 
for warm plasmas than for cold plasmas because the Vlasov equation has 
to be solved to second order in the perturbations. It follows that the 
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Lagrangian method should be even more advantageous for warm plasmas 
than for cold plasmas. No comprehensive development of the method has 
been presented so far for nonlinear wave interactions in warm plasmas. 

This is the principal aim of the present work. 

I.5.I The Microscopic Treatment 

Although a microscopic Lagrangian function for warm collision- 

62 

less plasmas was first developed by Low more than ten years ago, only 
a few authors have made use of it. There are certain obscurities in 
the original derivation. For example, Low showed that his Lagrangian 
gives the correct dynamical equations. In order to obtain small signal 
equations in correct form from the second order Lagrangian, however, he 
had to impose a consistency condition. Any solution derived from the 
expanded Lagrangian which would not satisfy this consistency condition 
was to be considered an incorrect solution. No discussion of the origin 
of the condition was given. 

Low' s Lagrangian has been used for several analyses of wave-wave 
interactions in warm plasmas: Suramlishvili obtained the coupling 

coefficients for three- and four-wave interactions involving Alfve^, 
helicon and ion acoustic waves Boyd and Turner^ considered the 

64 

interaction of two ordinary waves and a cyclotron harmonic wave. 

Dewar^ used Low's Lagrangian to show the important role of nonlinear 
Landau damping in the modulational instability. The emphasis of these 
contributions is rather on solving specific cases than providing general 
theories. Furthermore, they have not discussed Low's consistency condition 
for the second order Lagrangian, nor for those for orders higher than 
second. One may therefore question whether any of their solutions are 
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incorrect . 


This point has been cleared up by Galloway and Kim, 


who 


66 


have shown that the consistency condition is not necessary. 

1.3.2 The MHD Approximation 

So far, most applications of the Lagrangian approach to 
plasmas described by the hydrodynamic equations utilize the MHD approxi- 
mation, i.e., the plasmas are able to maintain a local isotropic velocity 
distribution by virtue of an assumed short mean free path for ion- 
neutral collisions. This model is used mainly to study low frequency 
phenomena that are usually associated with ion motions; the electrons 
serve primarily to maintain electrical neutrality. Higher frequency 
phenomena, such as plasma oscillations, are neglected. An immediate 
consequence of the restriction to low-frequency phenomena is the simpli- 
fication that the displacement current term in Maxwell's equations can 
be dropped. The current due to free electrons is supposed to obey 
Ohm's law, because of the high electron-ion collision frequency. 

The Lagrangian appropriate to the MHD approximation gives only the 
second moment equation (the momentum transfer equation); other moment 
equations, Ohm's law, and Maxwell's equations, are used as subsidiary 

frj 

conditions. This Lagrangian has been used by Eckart who showed that 

the equations of motion for both incompressible and compressible fluids 

68 

could be derived from variational principles. Katz derived the 
equation for a one-fluid, inviscid, perfectly compressible plasma in an 
electromagnetic field. Newcomb, ^ on the other hand, treated a perfectly 
conducting fluid governed by the conventional hydrodynamic equations, 
for which the displacement and convection current are negligible under 
most conditions of interest. Using the minimum potential energy 
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principle, he studied the stability of steady flow for isotropic 
plasmas in toroidal and poloidal fields. Similar systems consisting 

of fluid, vacuum and perfectly conducting solid walls were considered 

70 

by Lundgren. Application of this Lagrangian has been extended 

42 

recently to nonlinear problems by Dewar, who has studied the nonlinear 
interactions of hydromagnetic waves. 

1 . 4 Aim of the Research 

The foregoing review will serve to demonstrate the need to develop 
a unified and general theory of the Lagrangian method for nonlinear 
wave interactions in warm plasmas. This is the main purpose of our 
work, and theory will be presented for both the microscopic treatment 
and the hydrodynamic approximation to it . 

1.4.1 The Microscopic Treatment 

In what follows, we shall discuss first the microscopic 
treatment. This treatment is parallel to the averaged Lagrangian method 
developed in cold plasmas. A rigorous theory is presented to show how 
the method has to be modified in phase space. Results of some illustrative 
examples are compared to those previously obtained by other methods. 

In Chapter II, a Lagrangian for a warm plasma is formulated based on 
the work by Galloway and Kim.^ The origin of Low' s consistency 
condition is pointed out, and the lack of necessity for it is proved. 

This Lagrangian is used in Chapter III for the linear theory to obtain 
the equation of the particle trajectories and the dispersion relations. 

By proving the equivalence of the conservation of action and the conser- 

43 

vation of energy, the equivalence of the procedures used by Dougherty 

33 39 

and Dysthe, and by Galloway, for obtaining the coupling coefficient 



is shown. In Chapter IV, nonlinear wave-wave interaction is considered. 

The averaged Lagrangian method is applied to obtain the action transfer 
equation and the coupled mode equations. As some illustrative examples, 
the coupled mode equations are obtained in explicit forms for all 
possible interactions among waves propagating nearly parallel, and among 
those propagating exactly perpendicular, to the static magnetic field. 

These examples serve well to demonstrate the advantages of the method. 

1.4.2 The Hydrodynamic Approximation 

A plasma is well defined if the three velocity distribution 
functions for electrons, ions, and neutral molecules are known. However, 
they are often very difficult to determine. One must then be satisfied 
with certain mean values describing less perfectly the state of a 
plasma in the hydrodynamic approximation. This approximation is especially 
appropriate for a plasma where particle velocities are not comparable 
to the phase velocity of waves. 

Chapter V is devoted to the Lagrangian appropriate to the hydro- 
dynamic approximation. It is formulated by two different methods. The 
first is by integrating the exact microscopic Lagrangian in velocity 
space. For the linear and nonlinear theory for small perturbations, the 
Lagrangian is expanded. The expansion, however, involves some mathema- 
tical difficulties. The other method is to obtain the expanded hydrodynamic 
Lagrangian directly from the expanded microscopic Lagrangian by integrating 
in velocity space. In this way the difficulties involved in the expansion 
of the hydrodynamic Lagrangian can be avoided. As applications, inter- 
actions among two electron plasma waves and one ion acoustic wave, and 
among one electron plasma wave, one ordinary wave, and one ion acoustic 
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are considered by the first method, and that among two electron plasma 
waves and one Alfv^n wave by the second method. 

The Lagrangian to be described is more general than that of the 
MHD approximation mentioned in Section 1.3.2: first, it includes both 

the moment equation and Maxwell's equations; second, while only the 
ion motion is considered in the MHD approximation, our Lagrangian 
includes both electron and ion motions so as to account for both electron 
and ion waves; third, the MHD Lagrangian can be obtained from ours as 
the special case of spherical adiabatic compression. 
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Chapter II 


A LAGRANGIAN FORMULATION FOR WARM PLASMAS 

The complete set of equations for describing a warm collisionless 
plasma consists of the equation of motion for a charged particle, 

Maxwell’s equations, and the Vlasov equation. The dispersion properties 
of plasma waves are determined by combining the expressions for current 
and charge densities obtained from the Vlasov equation with Maxwell's 
equations . 

An alternative method for describing plasma waves is to use 
Hamilton’s principle. This can be done by formulating the Lagrangian 
in phase space to include the equation of motion, and Maxwell's equations. 
In the course of carrying out the variational procedure, a conservation 
law in phase space (Liouville's theorem) must be imposed on the evolu- 
tion of the electron distribution. For the collisionless plasma, this 
is the Vlasov equation. 

In Section 2.1, the exact expression for the Lagrangian will be 
formulated rigorously. It will be proved that this is the desired form 
by showing that the Euler equations give the equation of motion and 
Maxwell's equations. In Section 2.2, we shall deal with the Lagrangian 
for a slightly perturbed system. The required approximate form is 
obtained by expanding the Lagrangian about the equilibrium state. In 
Section 2.1, Lagrangian and Eulerian coordinate systems are adopted 
simultaneously for the particle part and the field part of the 
Lagrangian, respectively. In order to overcome the difficulties arising 
from using mixed coordinate systems, a field-like Lagrangian function 

37 

will be derived in Section 2.2, by making use of displacement vector 

71 

(sometimes called polarization vector ) coordinates. In Section 2.3> 
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it is shown that the first order Euler equations, obtained from the 
second order Lagrangian function, are the desired first order equation 
of motion and Maxwell's equations. Comments on Low's procedure are 
made in Section 2.4. 

2.1 The Field-Particle Lagrangian 

The Lagrangian for a system of electrons in an electromagnetic 
field consists of three contributions: that of the field, that of the 
electrons, and that of the interaction between the electrons and the 
field. In obtaining the last two, a Lagrangian coordinate system has 
to be adopted, while an Eulerian coordinate system is used for the 
first. We will first consider the last two parts separately. The 
total Lagrangian will then be obtained by adding the field contribu- 
tion to the electron contribution. 

The contribution of a single electron at > at time t > 

including its interaction effects with the field, will be the classical 

36 

expression for the charged-particle Lagrangian, 

= I mv &( t ) + < *[ < P( ( 2^, t ) “ • £(,*,* *)] t (2- 1 ) 

where q is the charge of an electron; cp and A are the scalar and 
vector potentials, and the subscript a indicates a specific electron. 
Equation (2.1) holds everywhere along its trajectory. 

Now consider an infinitesimal box, cell cr, in phase space. It is 
considered to contain many electrons in the vicinity of electron cr, 
and defined to move in time such that no electrons flow in or out, as 
shown in Figure 2.1. This conservation law implies that the phase 
volume and the electron density of cell a remain constant along its 
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Figure 2.1 Conservation in phase space [the number of electrons, n^ , 
inside cell cr is constant along its trajectory] . 
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trajectory. Since the electromagnetic fields can be considered uni- 
form over the cell, the contribution of cell a to the Lagrangian func- 
tion may be written as n^X^, where n^ represents the number of 
electrons in the cell, normalized to the total number of electrons in 
the system. The total contribution of the electrons in the system 

can be found by summing / over all a. The discrete summation of 

a 

cells can be transformed to integration by introducing continuous variables 

for c T, so that ^Tn^ ]d 3 r ' jd 3 v 'F(r',v' ) / , where F is a con- 

a J J 

tinuous function. The continuous set of variables (^r ,^v ) are known 
as Lagrangian coordinates, and are usually chosen to indicate the initial 

72 

position in phase space of cell a. 

F 

Adding the contribution of the field, X , to that of the electrons, 

P 

X , we obtain the total Lagrangian, 

L = Jd 3 ^ +J*d 3 r£ F , (2.2) 

P L 

where Z and Z are defined as 


Z = - mv 2 (r',v',t) + q[cp - v(r',v',t) • A] , 


1 

2 


Z 


F 


1 

2 



— (V X A 



(2.3) 


'Here, all quantities for which the independent variables are not ex- 
plicitly shown in X are to be considered as functions of ^r ,^v 
and t, through their dependence on the Eulerian variables r^, e.g., 
9 = cp[r(r , ,v , ,t),t] . 

The quantities to be independently varied are the position of an 
electron^ r, and the fields, A and all of which satisfy the 
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usual restrictions on their variation, 


6<r(t 1 ) = 6il»(t 2 ) = 0 , 


(2.4) 


where and t g are the initial and final times. Before taking 

variations, it should be noted that the conservation law in phase 
space requires 


dF 

dt 


= 0 . 


(2.5) 


Some authors have modified the problem by introducing a Lagrange 

multiplier for this subsidiary condition, but this is not essential. 

By letting F(r',v') = f(r,v,t), and changing from the variable 
set ( r 7 , v 7 ) to ( r, v) , we see that Equation (2.5) can be written as 


df _ df 
dt dt ~ 


F 

Vf +~ 
m 


V f = 
v 


(2.6) 


which is the Vlasov equation. 

With the help of Equations (2.4) and (2.6), the Euler-Lagrange 

equations from r, cp, and A variations of Equation (22), respectively, 
73 

yield J 


dv, 

m = -q(E + v X B) , 
at 


(2.7) 


where 



dA 

E = -Vcp - ^ , B = V X A . 


( 2 . 8 ) 

(2.9) 


( 2 . 10 ) 


Equation (2.7) describes the motion of an electron in the electromagnetic 
field. Equations (2.8) - (2.10) constitute the Maxwell equations. 
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2.2 Approximation to the Lagrangian 


In setting up a variational principle for the trajectories of an 
assembly of electrons, it is quite natural to use Lagrangian coordinates. 

As one may appreciate from Section 2.1, however, this would involve some 
difficulties in integrating by parts, since the differential operators 
are in Eulerian coordinates. Furthermore, the Eulerian coordinate system 
is overwhelmingly popular in describing plasmas. It is therefore con- 
venient to choose Eulerian coordinates for the independent variables. 
Fortunately, the conflicting requirements can be satisfied by the intro- 

37 

duction of a 'displacement vector'. 

Figure 2.2 shows a cell in phase space at (r,v), at an instant of 

time t, following a trajectory in a perturbed field. It would be at 

(r , v ) at the same instant of time if the perturbation were absent. Let 
/ "'-0 

( I, Jp be the displacement of the cell due to the perturbation, then 

Let (r , v ) be fixed for all time. The electrons occupying the cell at 
(rQ,Vo) in equilibrium at t are flowing out of the cell as time passes, 
and another group of electrons occupies the cell at a later time. The 
cell at (r,v) in the perturbed state, which corresponds to the equili- 
brium cell at (r , v ), is then located at a different place from the 
one at t. However, the new location (r,v) still satifies Equation (2.1l). 

From the definition of two corresponding cells, the number of elec- 
trons in the equilibrium cell is always the same as that of the per- 
turbed cell. Furthermore, the number of electrons remains constant for 
all time. It therefore follows that 

f (^t) d3rd3y = f o ( Io^ ) d V v o ' (2 - 12) 
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where f is the electron distribution function in the equilibrium 
state. But the conservation of density along the perturbed trajectory 
requires 

F(r',v') = i(r,v,t) , d 3 r'd 3 v' = d 3 r d 3 v . ( 2 . 13 ) 

Separating the fields into equilibrium and time varying quantities as 

<P . #0 + “fl ' £o + l ’ < 2 - l4 > 

and substituting Equations (2.1l) - (2.l4) into Equation (2.2) gives 
the Lagrangian 

L = J d \pVo<iW * P + J d3r ** ’ f 2 - 1 ?) 

where we have 

= t m (^o + 1 )2 + q f § 0 + ' C*0 + 8 • (&Q + * )] ' (2 - 16) 


£ F = 




$ 

0 0 


St 


Vi 


+ 


Sa\2 
St ) 


2|i, 


( V 0 X ^ + X £)‘ 


0 


~-0 


(2.17) 


The field quantities in Equation (2.16) should be evaluated at 
(r + v n + ?)> and the subscript 0 on the differential operators 
represents the (r , v ) - coordinate system. In Appendix A it is shown 
that 

5 = D? , (2.18) 

' r-< 


where the operator D is defined by 



^0 


7 o + ^0 


vO 


(2.19) 


The Lagrangian in Equation (2.15) is now a 'field-like* one, and the 
independent variables in Equation (2.15) can be considered as Eulerian 


coordinates . 
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If the perturbation is small, then the expansion of the Lagrangian 
in Equation (2.15) generates a series in which L r , the Lagrangian of 


n th order in the perturbation, is given by 


L n = J d3r / d3v f ofe)^ + J d3r 4 • 


( 2 . 20 ) 


Here the subscript 0 on r and v has been dropped, and, unless spe- 


dfied, the subscripts 0 and 1 will be dropped from variables hereafter 
up to Chapter V. The first few terms of and are as follows: 

Zero Order: For the unperturbed terms we have 


= i mv 2 + q($ - v • A) , 


T) 


= T ( ™ )2 ‘ 2^ (7X V 2 • 


( 2 . 21 ) 


First Order: We have the separate expressions. 


= mv • D£ + q[<p + (J£ • ^)§ “ X ’ (,£ " - ' t ” Z ' 0) > 


h m e o 74 ' 


+ It) - jr: < v x ’ (vxa) 


( 2 . 22 ) 


It is obvious that 


L 1 = 0 , 


(2.23) 


since K dt has been made stationary by choice of the zeroth order 
solution . 

Second Order : In this case, the separate expressions are 

4 = \ m (D|) 2+ q[( 5 - *)<P+| (|§:W)f]-qv. [(§'V)a+| (S|:W) 2 A] 


-* D i • t* + (i • V )A] > 4 = + I?) - 2^ (vx ~ )2? (2> 
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where the double dot product is such that ab:cd = a . (b • c)d. Sub- 
stituting them in Equation (2.20), and making use of the relations, 

e = ,-VCp - , b = V X a , (2.25) 

gives, after some manipulation, 

1 2 ’ J d3r j d3vf oii m(D?)2 - “I • <® + x x £> + 1 " D I ■ I x ® 

- | q i • (I • V)(E + v X B)J + ijd 3 r (e 0 e 2 - i b £ ) . (2.26) 

Third Order: The field contribution vanishes at this order, leaving 
4 | q[( 55: W)cp + i (S5S-W7)*] - \ qv • [(5?:W)a + \ (K5:W7)AJ 

-q D? • [(? • V)a + i (??:W)A] , Z = 0 ' (2-27) 

e± r^r^t r^> J 

where the triple dot product is defined such that abc:de £ = a • [ b- (c -d)e]f . 
Substitution in Equation (2.20) gives 

1 3 = q f d3r j" d3v f o{l (?£ :W ) Cp+ ^ (£||:w*)*-v .[| (H?-^)£I 

- D? • [(5 • V)a + \ (?5:W)A]{ . (2.28) 

2 . 3 First Order Equations 

Before proceeding further, we must demonstrate that variation of 

our second order Lagrangian gives the small signal equations correctly. 

Variation of I-: Taking the variation of L [Equation (2.26)] with 

'y 2 

respect to ^ yields 

mD 2 ? = f + (F • V)F - q X B , (2.29) 

where F and f are given by 
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(2.30) 


F = -q(E + v X B) , f = -q(e + v X b) . 

Now, expansion of Equation (2.7) about (£ q >^ q ) i n equilibrium yields 

m 3t <£o + 8 - " q [£ + (i • 7 )£ + £ + + i) x f£-K£- £)£ + £}]- < 2 -3 1 ) 

Extracting the first order quantities from Equation (2.3^)> an d dropping 
the subscript 0 from v^, gives Equation (2.29)> and thus confirms that 
it is the correct small signal equation of motion. 

Variation of <p: Carrying out this variation leads to the expression 

e 0 V ‘ ~ = q | d3v7 * ( f o£) > (2.32) 

which we must show to be the first order Poisson equation. To do this, 
it is sufficient to show that the first order charge density is 


P 


‘I 


d 3 vV 


( f o£) 


(2.33) 


From the definition of a cell, described in Section 2.2, it follows 
that the number of electrons in the cell at (r,v) in the perturbed field 
is the same as that in the cell at (r^v^) in the equilibrium state 


f(r,v,t)d 3vd3r = f 0 (^,^)d 3 v o d 3 r o , 


(2.3*0 


where d 3 vd 3 r, the volume of the cell at (r,v) (see Figure 2.3)> is 
related to d 3 v Q d 3 r Q , the volume of the cell at (r^Vg)* by 


AA . a VA„ , j = , 

00 


(2.35) 


where 3 is the Jacobian. Substitution of Equation (2.1l) for the 
Jacobian reduces Equation (2.35) to the first order equation 


AA = (1 + V Q • _|)(I + V v0 • i)d 3 v Q d : 


0 


(2.36) 
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Figure 2.3 Phase space elements [the volumes of the shaded ele- 

3 3 

ments are both d r d v . The volume of the cell in 
the perturbed system is d^rd^v] . 
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From Equation 
cell at (r, v), 



(2.34), the number of electrons contained in the perturbed 

3 3 

of volume d r^d v^, may be written to first order as 


d3v n d3r n = - 7 n • E - v • ?)dVd J r . (2.37) 


3., „3, 


0 




vO 


0 


On the other hand, the number of particles in the cell at (r,v), of 

3 3 

volume d r^d v in equilibrium, is to first order 


f ote )d3 V 3r o = [f 0 ( -Jo^ ) + ( £‘W^^o )+ ( i' 7 vo )f o ( ^^ )]d3 V 3r o • 

(2.38) 


Subtracting Equation ( 2 . 38 ) from Equation (2.37) gives the first order 
distribution function 

‘itez) - - 7 o • <VS> - v vo • (f o9 • < 2 - 39) 


Integrating f T over v_ yields 
1 

p = qJd 3 v 0 V 0 • (IqD • (2. hO) 

Dropping the subscript 0 from ( r • , v.) in Equation (2.40) gives Equation 
(2.33)- Equation (2.32) is, therefore, the correct first order Poisson 
equation . 

Variation of a: Finally, variation of Equation (2.26) with respect 

to a, and use of the conservation law for the equilibrium density, 


v • Vf + a • V f = 0 , 
~ 0 ~ v 0 


(2.41) 


yields 


jr v » £ = «o st - - X’ ■ • ( 2 - h2 1 


"0 

It only remains to show that the first order current density is 


A= ->pV £ ol - Zo v o • < f oS>] > 


( 2 . 43 ) 
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to confirm that Equation (2.42) is one of Maxwell's equations. 


By following an argument similar to that used above, we may write 


v f(r,v,t)d 3 vd 3 r = < < 2 ' 44 > 


j3. 


for the perturbed cell at (r,v). Substitution from Equation (2.36) gives 


v l(r,v,t)d 3 v 0 d 3 r 0 « l L:] + .C)t 0 fi3,;i 3 )d 3 v 0 d 3 r 0 


-I + V v0 • a d3v d d \ • (2- 4 5) 


0 0 


The term similar to Equation (2.44) for the cell at (r,v) in the equili- 
brium state may be written as 


Z f o^ d3 V 3r o = '!Ui + I) f o d3v o d3r o 


+ *<£ • 7 o + S ' ■ f 2 - 46 ’ 


The first order current density is 

j = q(d 3 v v[f ( r, v,t) - f (r,v)] . (2.47) 

Substitution of Equations (2.45) and (2.46) reduces Equation (2.47) to 
Equation (2.43), as required. 


2.4 Comments on Low's Procedure 

Low (1958) expanded the Lagrangian about the equilibrium state 
according to the scheme illustrated in Figure 2.4. Because of the ini- 
tial conditions 

? = S = 0 , (2.48) 

r*y 

O O 

the initial electron density, F(r , ,v / ), and volume, d~r d'V ; , in the 
perturbed system can be considered as those at t = 0 in equilibrium, 
as shown in Figure 2.4. The conservation law along the unperturbed 
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Figure 2.b Low's model for the expansion of the Lagrangian. 
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trajectory then gives 


F(r',v') = f(r,y,t) = = d\d 3 v n , (2.49) 


0 0 


and it is this equation which Low used for the expansion of the Lagrangian 
Equation (2.49) requires unit Jacobian, which use of Equation (2.1l) shows 
is equivalent to 


• E + V • § = 0 . 

V ~ 


(2.50) 


This is the so-called "consistency condition", used by Low in proving 
the first order Euler-Lagrange equations [Equations (2.32) and (2.42)] to 
be the first order Maxwell equations. Low argued that any solution from 
his Lagrangian which does not satisfy Equation (2. 50) must be thrown out 
as an improper one. 

If the consistency condition is accepted, the misleading impression 
is gained that the Lagrangian would be invalid in a system with different 
initial conditions, or in a steady state which is not affected by any 
initial conditions. Furthermore, a higher order "consistency condition" 
would presumably be required in order to check the validity of solutions 
to nonlinear problems, i.e. corresponding to the inclusion of and 

higher order terms in the series approximation to £. 

On the other hand, our transformation of coordinates via Equations 
(2.12) and (2.13) does not require any specification of initial conditions 
Neither do we need any restriction in proving Equations (2.32) and (2.42) 
to be Maxwell's equations. Low's Lagrangian [Equation (2.2)] and the 
supplementary conservation laws [Equation (2.5)] are in fact the complete 
set of equations to describe a warm, collisionless plasma with arbitrary 
initial conditions. All solutions from his Lagrangian are therefore 
proper without invoking his consistency condition. 
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Chapter III 


MICROSCOPIC TREATMENT 
I . LINEAR THEORY 

To describe the small-signal behavior of warm plasmas, we may 
use the first order Euler-Lagrange equations [Equations (2.29)> (2.32) 
and (2.42)] obtained from the Lagrangian second order in perturbation 
in Chapter II. Fourier transformation of those equations should give 
the dispersion relation. An alternative procedure is to average the 
second order component of the Lagrangian density with the perturbing 
variables expanded in a Fourier series. Variations of the averaged 
Lagrangian with respect to the Fourier amplitudes of the variables give 
the Fourier transformed form of the first order equations. The advantage 
of the averaged Lagrangian method is that it gives an additional equation 
which cannot be derived in Chapter II. It is the equation for conserva- 
tion of action obtained from the variation with respect to phase. 

The above procedure is presented in Section 3-1 » small signal equa- 
tions and the equation for conservation of action are obtained for an 
infinite, collisionless, warm plasma in a uniform static magnetic field. 

It is also shown that the Euler-Lagrange equations obtained by variation 
with respect to the amplitudes are the first order equations obtained in 
Chapter II. In Section 3-2, the equation for the conservation of energy 
is derived from the averaged Hamiltonian, and is shown to be equivalent 
to conservation of action. The particle trajectory equation is solved, 
and wave dispersion relations are obtained in Section 3-3 for all possible 
magnetoplasma waves propagating at an oblique angle to the static magne- 
tic field. Results of Section 3.3 are given for the limiting cases of 
waves propagating nearly parallel and exactly perpendicular to the magne- . 
tic field in Sections 3.4 and 3 . 5 . 
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3.1 Euler-Lagrange Equations from the Averaged Lagrangian 


3.1.1 Procedure 

The perturbing dynamical variables can always be expressed 
as a sum of propagating waves 

?(r,v,t) = E'^v) exp[i0 k (r,t)] + §^ k (v) exp[ i0_ fc ( r, t ) ] } , 
k 

a(r,t) = exp[i9 k (r,t)] + a^ k exp[i0_ k ( r,t)} , (3-l) 

k 

where a prime indicates that the summation extends over only half of 
k-space, and the 0 's are integer multiples of a phase angle 9. 

rv< K 

Then variables ST and a, satisfy 
~k ~k 

t) exp(-ik0)d0 , 

~k = f ex P( -ik0 ) d0 • (3.2) 

T) 

An equation similar to the second one in Equation ( 3 . 1 ) can P e written 

for the scalar potential, cp. It is apparent then from the Fourier analysis, 

that 



4k = 4c ’ ~k " > "-k “ "k 


e , = -0,. • 


(3.3) 


The angular-frequency, u^, and the wave vector, k, are given by 


a6 k 

^ = at" > 


k = - 


^k 

ar • 


( 3 .^) 


The and k are com P lex i n general. However, their imaginary 

parts can be assumed to be very small compared with their real parts for 
stable or nearly stable waves, for which the linear theory is valid. The 
imaginary parts may therefore be neglected. 

With the definition of Equation (3.1), the single particle Lagrangian 
second order in perturbation may be written in the form 





(3-5) 




(<V k i>5 k >4 



) 


1 


where the bar over indicates that averages are taken in space and 

j* . . . 

time. The Euler-Lagrange equation obtained from the variation 01 

Equation (3-5) yields 


s/( 2 ) 

k 




= 0 , 


(3.6) 


Solving Equation (3.6) gives the spatial displacement vector, j^, in 


terms of field quantities. 

The averaged Lagrangian density component, Z , is also a function 
of # However, since can be eliminated with the help of Equation 

(3.6), it can be written as 


T = s' £ 

2 k k 


( 2 ) 



Z 





i 




(3.7) 


j* . 

The Euler-Lagrange equations obtained from the a k and 0 k variations 
of Equation (3.7) yield, respectively, 


a 

St 



(3.8) 


(3-9) 


The derivative B/Bt (or 3/3r) denotes that the differentiation is taken 
with r (or t) fixed. It is, however, considered as 'total' in the 
sense that the dependence of on t (or r) through or k is 

allowed for. The nontrivial solution of Equation (3.8) gives the 
eigenvalues determined by 

det = 0 , (3.10) 

K 
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(3.H) 


Its eigenvectors reduce Equation (3.7) to 


. 0 . 


Equation (3.9) may be written as 


S_ 

dt 


(\) 


a 

¥7 





(3.12) 


where 



(3-13) 


Equations (3.6) and (3.8) are the electron equation of motion and 
Maxwell's equations, respectively. Equation (3.10) is the small signal 
dispersion relation. Additional information is provided by Equation 

(3.12) , which could not be obtained in Chapter II. It can be shown that 

48 

N is the action variable for a system with periodic motion. Equation 
k 

(3.12) is therefore to be interpreted as the conservation of action. 

3.1.2 Explicit Expressions for the Euler-Lagrange Equations 

Explicit expressions for the single particle equation of 
motion [Equation (3-6)], Maxwell's equations [Equation (3-8)], and the 
equation for the conservation of action [Equation (3.12)] are required. 
However, that for the action conservation equation is not given below. 

It will be shown instead in Section 3-2 that the conservation of action 
is equivalent to the conservation of energy and momentum. Since the 
energy equation is more familiar to plasma physicists, the explicit ex- 
pression for it is derived in Section 3-2. The explicit expressions for 
Equations (3-6) and (3-8) can be found by making use of the second order 
component of the Lagrangian function obtained in Chapter II. 
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Let a plasma in a uniform static magnetic field oriented in the 
z-direction be subject to a small perturbation. The zero order quantities 
may then be assumed to be homogeneous . With the gauge chosen such that 
$ = cp = 0, averaging the Lagrangian functions in Equations (2.24) and (2.26) 
gives the single particle Lagrangian and the Lagrangian density, respec- 
tively, as 


£ k (2) = i m |D k,Al 2+ I q(D k,oik) *4 x £ +q 4-(^ + ^ x 4) + c - c - (3-^) 


4 2) = " I q j* d3y f o % * <4 + z x V + \ 


" 2^ 0 ^ + c,c *> (3-15) 


where c.c. indicates complex conjugate and 


e. = - ioi , a. , 
-k k— k ’ 


4 = x 4 > 


( 3 . 16 ) 


D S <n 
k,p c 


Ip + 1 - Vi * <VVz + '“c « 


(m- = 0, ± 1 ,± 2 , ...). (3.17) 


In Equation (3*17)> 00 , is the electron cyclotron frequency; ({) and (j)^ 
are the azimuthal angles of v and k, respectively, and the subscript x 

/VI /VI 

indicates a component perpendicular to the static magnetic field (see 
Figure 3-1) • It may be seen immediately that D is the Fourier trans- 

K y U 

form of the RHS of Equation ( 2 .I 9 ). We note that we may set 4=0 without 

K 

loss of generality. However, in nonlinear problems involving interactions 
among many waves, all of the waves do not necessarily propagate collinearly, 
i.e., the <f> k 's of all waves are not zero. In preparation for the 
following chapter, we shall assume <j) k £ 0 here. After substitution of 
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X 



Figure 3*1 Coordinate system for particle velocity 
and wave vector. 
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of Equations (3-l4) and (3.15)> Equations (3-6) and (3-8) become 


D k,0 4k + °k,0 4k X 


CD = - — 
~c m 


Yi-£=L ♦ *-<» 

“kT \ ~ 


a,) 


~-k 


(3.18) 


B x <B x 4k) + • 4 k = 0 * 


(3-19) 


where ft ( = kc/tn, ) is the refractive index vector and cd is a vector 
with its magnitude co c directed along the z-axis. The equivalent plasma 
permittivity tensor, k^, and the current density, j^, are given by 


3. ' 1 + - - 1< >j d3vl o [ “k t ik + £ * <i x • (3 ' 


£ o“kSk 

Equations (3.18) and (3.I9) are the equation of motion of an electron and 
the wave equation, respectively. It should be noted that these equations 
are Fourier transforms of Equations (2.29) and (2.42). 
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3.2 Energy and Momentum Conservation 

To each dynamical variable, Y^, corresponds a canonical momentum, 

it , defined by 

M* 




- 9(9Y /at) • 

H* 


(3.2i) 


However, no canonical momentum exists for the Fourier amplitude of the 

i ( 2) / \ 

dynamical variables a^ or a^ , because ' in Equation (3-7) does 

not involve a or a^ . Only the phase, 9 , has a canonical momentum. 

K K K 

This is the action variable, N [defined in Equation (3-13)] • 

t- XV 

According to Hamilton’s theory, and with the use of Equation (3.11), 
the generalized energy and momentum densities averaged in time and 
space are obtained as 


*2 = l' 4 2) - 


K ( 2) = <\"k 


(3-22) 
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(3-23) 


& = [' si 2 ’ - si 2) = i \ 


( 2 ) 

If we differentiate ' and substitute Equation ( 3 . 13 ), we obtain 


d<cf 2 ) /3o5 3k\ /den Sk \ 

__iS — = n ( — — — - v • — 2 ) _J£ — = N ( — - _ v • _H J 

ot k \ot ~k 3t/ ’ or” k ySr^ — k Hrj’ 

By combining this with Equation (3-12), and the ray equation obtained 
from Equation (3.4) 

3k 3o5 


5T = " 



> 


(3-24) 


(3.25) 


we may express the total differentiation of 
to time as 



and 



with respect 


|- (05 N ) 

3t 1 k k' 


3 

17 




05 


N ) = 
k' 


( 3 - 26 ) 


It t k iV + tr • (Sk k i V - 0 • 

As in Equation (3-9); the partial derivatives in Equations (3.24) - (3.27) 
merely indicate that the differential is taken with other coordinates 
fixed. They are, however, considered as ’total’ in Hamilton’s theory. 
Equations (3-26) and (3.27) are the equations for conservation of energy 
and momentum, respectively. It is important to note that the use of 
Equations ( 3 . 11), (3*24) and (3.25) reduces the conservation of energy 
[Equation ( 3 . 26 )], or momentum [Equation (3.27)], to the conservation of 
action [Equation (3.12)]. 

From Equations (3-13), (3-15) and ( 3 .I 9 ), explicit expressions may 
be obtained for the energy and power flow: 


05 N 

k k 


_1 


b. 1 2 + - 


0 


~-k 


2 ~-k 


’0 $05 


K 


K.) • 

~k' 


•^k 


+ c.c. 


(3.28) 
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1 


(3.29) 


“k "i. = ^ 4 x i£ - ! S ( 4 


£ 0 “k ~k 


-%) + c - c - 


The first term on the RHS of Equation (3.28) is the magnetic energy in 
free space. The second contains the electric energy in free space, the 
particle energy, and the particle-field interaction energy. The first 
term on the RHS of Equation (3.29) is the Poynting vector, and the second 
is the power flow due to particles flowing out of a volume coherently 
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with the wave. 


3.3 Solutions of the Euler-Lagrange Equations 
3 . 3 .I Perturbed Electron Trajectory 

A difficulty in solving Equation (3-18) is that the x and y 
components of the LHS are coupled. It may easily be seen, however, that 
it can be diagonalized by a unitary matrix 



i 0 
-i 0 
0 v/2" 


(3.30) 


New characteristic vectors for are then 


<£> & 4 )> where 


5 k = zo ij ! k> j = *.?,*)• 

J 


(3.3i) 


Since the y component is equal in magnitude to the x component, and 
o R 

lags it by 90 , 5 " is a vector rotating in a right-handed sense about 
K 

the z-axis. Similarly is rotating in a left-handed sense. The 

unitary transformation of Equation (3.^8) gives 


D, • E. = M, • e , 
~k '■5k ~k ~k ' 


(3.32) 
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where we have 


4 = 


D D 
k,-l k,0 

0 

0 


°k,l D k,0 

0 


0 
0 

Vo 


(3-33) 


m = 

~k 


m RR 

M k 






, ZR ZL 

Lm, m. 


,RZ 


m: 


M, ZZ J 


(3-34) 


The elements of M, are given by 
~k 


RR 


M 


= «“■* m ,s(i« . 




k v 
i l 

2co 


exp L i (<|)-<1> k ) ] | , 


M. 


ZZ 


- - 2 D - <=•• <K>1 . 


k k m 


exp [i((j)+<j) )] , 


,.ZR ZL* 

\ = \ 


k v 

q z x , 

— — exp( -i( 

m , 

V2 * 

k 


m; 


RZ 


= yt z * - - a 


a V* 

m 42 a> 


exp(i<j> .) . 


(3-35) 


After some manipulation, the solution to Equation (3.32) can be 
written explicitly as 


4 = 



4 > 


(3.36) 
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where the elements of the matrix G k are 


p o,-x F 0 (1 


( k v 

G zz = a P n — - 

k m) 0,0 2^ 


~[ pTl -l -l ex P[ i ( < M ) k ) ] + P l,l e x P[ i (4>-4> k ) ] 


k v 


= fST = * ^ , expLiCMj] , 


,RL 

\ ~ u -k m 203, ‘0,1 

k ’ 


ZR 

V 

-r 

= SL 

m 

k v 

Z 1 

n 

p -i,- 


LZ* 


k v 


RZ 

11 

3 

X Z 


G k 

= G -k 


P 

0,-1 


exp 


(- 4 ) , 


(3-37) 


n 


and P is defined by 


- \4 k > expjlWH.) + \ CH k )]j -(3.38) 

n=-oo k >n-*p k,n+v ' 


Here, X k is written for k^v^/oo,, and the following Bessel function 
identity has been introduced: 


exp (i^sintj)) = Y] J n U k ) ex P (in<|>) . ^ 3 ’ 3 

n=-oo 

The subscripts on the G* s in Equation (3.37) indicate the k or -k^ wave, 
as shown in Figure 3.2. A change in sign of the subscript k makes sign 
changes in k^ and k^, but the sign of the azimuth, <j> k , does not 

change . 

Equation (3.36) describes the perturbed electron trajectory. It 
may be seen that, if the rf electric field has a transverse variation 



X 
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?>t 


(k^ ^ 0), strong interactions occur between the rf field and the elec- 
tron when the Doppler-shifted frequency approaches a cycloton harmonic. 
3.3-2 Dispersion Relations 

X V z 

The unitary transformation of e = (e , e"\ e ) to 

K K K 

R L Z 

e fc = (e , e , e ) reduces Equation (3-19) 

K R K 


A • e = 0 , 
~ ~k J 


(3- 4o) 


where 


-n 2 -ji 2 /2 + « RR 
z x k 


-(71 2 /2) exp( 2i<|) k )+ ( n z 7l/\l2 ) exp( i<j> k )+ 


RZ -1 


-( Jl 2 /2) exp( -2i<j> k )+ k^ R -7l?-7l 2 /£t^ L 


z x k 


(Jl z T| i /72)exp(-i^ k )+K 


LZ 


(n z Tl i />/2)exp(-i<}) k ) + K k R (n z 7l x /\^)exp(i(|> k )+/< k lj -71* + < 


ZL n 2 . „ZZ 


(3.41) 


Here, 71 and 71 are the components of the refractive index vector 
z x 

parallel and perpendicular to the z-axis. , Using Equations (3.20) and 
( 3 . 36 )j and the relations 


J n<V - - r J „(V - f W - J „ + l< V > (3X2) 


2it 

I exp[i(m(|) + X k s in(j))] d(j) = 2 * J.„( x k ) , 


(3.43) 


the equivalent permittivity matrix can be obtained explicitly as 
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k - k H, 


y/_L HU. ( V/ , <Vi )2 l ^, \o k A JT \ 

+-\“k,n-lK,„ a^.n-l 2l \, n+ lJ " \„\ M1 » ” +1 / ’ 


,zz 


i _ n\ 2 E /r/v^\ 2 + ._ <y«> 2 1 A , 

W/ n\L\%n/ %,n+l %n-lj n / 


k rl k lr* 

k k 


/v 

vv 


exp 




^n+1 ^n-1 


(k v y 

V X x' 


%0 Vi 

^,11 ^^n+l 


J J 
n n +1 


2 %n+l %n-l n 
> ’ 


k 


= \ 


ZR* ZL LZ* 
= K = K 
-k -k 


- -(^) exp < 1, M£ 


‘VoVz 


«T - 


K +n “c )k z V x 


"n ^2 <\^ n ^2 q^ n 


J J ^ 
n n +1 j 


, (3.44) 


where is the electron plasma frequency; ( ) denotes the mean value over 

velocity space, and the argument of the Bessel functions is X . Unless 
otherwise indicated, the summations over n contain all negative and posi- 
tive integers. The small-signal dispersion relation is expressed by 
det £ = 0 [Equation (3.10)]. 

The solution of the perturbed trajectory of an electron [Equation 
(3.36)], the dispersion relation [Equation ( 3 . 40 ) ] , and the equations for 
the energy conservation [Equations ( 3 . 26 ), (3.28) and ( 3 . 29 )] are very 
complicated. Simpler forms can be obtained for the special cases of waves 
propagating nearly parallel or exactly perpendicular to the magnetic field. 
These two limits will be considered in the following sections. 
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3.4 Quasiparallel Propagation (k « k ) 

-L z 

3.4.1 Perturbed Electron Trajectory 


To first order in k^, Equation (3*37) reduces to 


,RR „LL* 
r k = G -k 


q 

"‘VVi 



k v 
x 1 



exp[ i(<j)-<|> k ) ] 


+ 


k v 

x x 


%-2 


exp[ — i ( <{) — 4> k ) ] 





4 


( k v 

ll + 11 

'VVcl-lv/' 

l 2 “k 

2 

L <V-1 J 


Lxp[-i(<H>. )] 


+ 


<\K +a >c )-(k z v z ) 2 
2 


ex P[ i ( <t>— <l> k ) ] 


y 


„RL ^LR* 
G, = G 
k -k 


q_ 

2m 


k v 

x x 


°V¥,tf¥,l 


exp[i(<jH-(j) k )] 


> 


ZR ZL* 

G, = G . 
k -k 


k v 

q 2 x 


exp(-i<j>) 



k v 
x x 




<0 



exp[i((j)-(J) k )] 


+ 


Vl (t \,-I + °‘k,-2 ) 



ex P [ -i ,l t i -'t i k ) 1 


> 


RZ LZ* 

G, = G 
k -k 


q 


k v 
x z 


2* m ^k 0 ^, 0^,-1 



(3.45) 


We shall make use of these expressions in Section 4.3. 


In the case of all waves propagating precisely parallel to the 


magnetic field (k^ = 0), Equation (3-36) reduces to 


s R = 




f = —3 

^ "Vk,i 



b 

-p ^2 

42 m 


+ 


i -±- ^ exp( i(J)) + -X_ 
' l2mm k,l "^,0 


(3.46) 
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The displacement of an electron in the xy-plane coincides in direction 

with the electric field of the right-hand or left-hand polarized wave, 

which interacts strongly with electrons traveling such that CD^-k^v^ ** 

± cd . The first two terms of the displacement in the z-direction are 
c 

due to the Lorentz force of the transverse magnetic field ( « v X b, ) . 

v ~L — k ' 

However, averaging over velocity space causes these terms to vanish. The 
average displacement in the z-direction is therefore only due to the 
longitudinal plasma (Langmuir) wave, which interacts strongly with elec- 
trons traveling with velocities close to the z-component of the phase 
velocity, i.e., oo^-k^v^ *=* 

3.4.2 Dispersion Relations 

To first order in k^, the elements of matrix A in Equation 
(3.4l) become 



.RZ ZR* a ZL .LZ* 

A = A, = A , = A , 
k k -k -k 




k V 2 

z 1 


'“k,o 


In the limit of purely parallel propagation (k^ = 0), the off-diagonal 
terms of A vanish. The nonzero diagonal terms are then the dispersion 

relations of the right-hand and left-hand polarized transverse waves, 
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and the longitudinal Langmuir wave. 
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3.5 Perpendicular Propagation (k = 0) 

Z 

3 . 5 .I Perturbed Electron Trajectory 

For this case, a rectangular coordinate system is more con- 
venient than the rotating coordinate system used in Sections 3-3 and 3.4. 

. R L Z. 

Putting k z = 0, and transforming the electric field (e k , e^, e k ) to 
(e|L e|~, e ), as indicated in Figure 3-3> reduces Equation (3-36) to 

K K K 


where 





(3.48) 



q 

m 




q 

m 



> 



»^ ± ( p o,-i e * p [‘‘(H)i + p o,i 


> 


G -MI --iS P n 

k m 0,1,-1 ’ 


<■? - - » * 


5 “ - G k' + 1 m { P 0,-l expC-iC^-^)] - P"^ exp[i(4.4 k )]J, 


_ k v 

of - a -±1 p" , , 

k m 1,-1 ’ 


_ k v 

of .-iSiir" 

k m 


<\ 0,1, -1 * 


(3.49) 


and, similar to Equation (3.38), P is defined as 

* y J ” 0,p,v 


s E ( ' l)n J ” (v) exp { i[n(H * )+ v^H)]{ 

n= '" (3.50) 



X 
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It should be noted that superscripts X and || indicate the components 
perpendicular and parallel to the wave vector, respectively. On the 
other hand, the subscripts x and || represent the components perpen- 
dicular and parallel to the static magnetic field. 

To see how each component of the rf field contributes to § in 

r* w* 

Equation (3.48), let us first consider the effect of the longitudinal 
component (ej^). The displacement of !jj[ due to ej| is related by 
jrll _ G ll II Jl tlie expression for which corresponds to that previously 

K K K 

75 II 

obtained by Crawford. This shows that a perturbation due to e" becomes 

very large when approaches (n ± 1 ) cq c , i.e., when the wave provides 

space harmonics of the electric field over the direction of the electron 

cyclotron motion. It may be seen that the contribution of ej^ to 

is due to an effect similar to the electron drift in a static magnetic 

field (v ~ e]^ x B/B^).^ The electric field component ej| makes no 

Z JL Z 

contribution to E^. The perpendicular components ( e^ and e k ) exert 

effects on ^ similar to those of ej^. However, they have perpendi- 
cular rf magnetic fields. Therefore, additional contributions are made 
by the Lorentz force (a v X b ). This can be seen by noting that the 
terms in Equation (3*49) giving these additional contributions all in- 
volve k^v^/oo^ or k^v^/o^ . With the help of Equation ( 3 . 16 ), k^v^e^/cu^, 

z Z Z x ’ X 

k v e, /oi and k v er/oi can be rewritten as v b , -v b, and -v b , 
xxkx xzkK x k x k z k 7 

respectively. 

3.5.2 Dispersion Relations 

For this case, the elements of matrix A are obtained from 

11 1 z 

Equation (3.4l), with the base vector (e", eT", e ), as 

K K K 
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= 1 


= 1 


2 Hr w 1 v ( v 
A W/ 
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<(V»> 2 > 


J J . 
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^^n-l^k^n+l 
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xil* ,2 




“¥“0 E 
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^An-l^n+l 


.xZ 


= A 


zx* .. a V U c k x £ 


= 1 


<v»> 


= A 


Z|| 


2 , 

^x 


“"^k n ^^n^^n-l 0 ^, n +1 

(v J 2 \ 

■> \ z n/ 




(3-51) 


The elements of the matrix A for perpendicular propagation have been 

/-N^ 

77 

obtained previously by Baldwin etal. Their results can easily be 
reduced to Equation ( 3 . 51 ) by making use of suitable Bessel function 
identities [Equation (3.42)]. 

For an electron velocity distribution function even in v^, the 
elements A^ Z , A Z-L , a||_ Z and A Z ^ all vanish, and e Z is decoupled from 

K K K K K 

II x Z 

e, and e, . The electric field e, is therefore a normal mode whose 
k k k 

ZZ 

dispersion relation is determined by A =0. In the cold plasma limit 

K 

this reduces to the well known dispersion relation of the ordinary mode. 

ZZ 2 2 2 

For a nonzero temperature plasma, however, A has poles at cd = n co ; 

K C 

n = 1 , 2,3 •••} which lead to propagation in the neighborhoods of the 

2 

cyclotron harmonics. For this reason, we shall refer to e, as the 

7 k 

ordinary cyclotron harmonic wave. The case where longitudinal waves 
(ej^) are almost normal modes has been considered in a great detail by 

77 

Baldwin etal . By this we mean that there exist solutions of the 
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dispersion relation other than ef; namely a||^ A XX - A^ 1 A X ^ =0. To 

K K K K K 

a good approximation these are the solutions of A^ = 0 and are 

JK. 

known as longitudinal cyclotron harmonic waves (or Bernstein modes). 

For the special case of an isotropic Maxwellian electron velocity 
distribution where 


/ » \ 3/2 

r / 2 ^ 2^ n 

m(v + v ) 
v j. z' 

\2jtkT ) 6XP 

2kT 


(3-52) 


one can readily obtain, upon substituting the integration identity 
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00 / 2 

f XJjj(X) exp/- dx = a exp(-a) ija) , (3-53) 

*0 


into Equation (3>5l)> the dispersion relations for the longitudinal and 
ordinary cyclotron harmonic waves, respectively in the forms 


1 = 


exp(-k 2 KT/mo> 2 ) ® n 2 I (k 2 KT/mn 2 ) 
t— x c i n x c 

2a v . 2 „ / 2 — E 


p p 

k KT/ma) 
l c 


^ CU k,n *\,-n 
n=l ’ ’ 


1 - 


n 2 I (k^T/mx?) 

n ' 1 C 


? ? 2 ? ” n vV T/I 

Ttf = ( — ) 1 + 2 a > exp( -k kT/wco) 

X \“c / C X C ^ %n%-n 

u n=l ’ ’ 


( 3 . 5*0 


Dispersion curves corresponding to the first and second expressions in 
Equation (3-5*0 have been computed by Crawford^ and Lee^, respectively, 
for a range of plasma parameters. The second term on the RHS of the 
dispersion relation for the ordinary cyclotron harmonic wave is negli- 
gible except very near to cyclotron harmonics, i.e., at nearly all 
frequencies the dispersion relation is essentially that of the ordi- 
nary wave in a cold plasma. 
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Chapter IV 


MICROSCOPIC TREATMENT 
II. NONLINEAR WAVE-WAVE INTERACTIONS 

The linear theory of Chapter III assumes that an arbitrary perturbation 
can be expressed as a superposition of non-interacting eigenmodes. A 
characteristic result of nonlinearity is the occurrence of coupling among 
these modes. If the coupling is weak, then arbitrary perturbations may 
still be expressed approximately as a superposition of eigenmodes, but 
whose amplitudes are now assumed to vary slowly in time and space. By 
'slowly' we mean that the scale of the variation due to the wave-wave 
interaction is much longer in time and space than the period and wave- 
length of the disturbance. We shall assume that the ratio of this scale 
is comparable to the ratio of the amplitude of a perturbation to that of 
the zero order quantity, and denote them by a small parameter e. We 
shall also assume that interaction occurs coherently, i.e., we shall not 
treat interaction in the random phase approximation described in Section 
1 . 2 . We shall consider all possible types of three-wave interactions 
in a warm, collisionless magnetoplasma, with an immobile ion background. 

The procedure to be followed is similar to that used by Dysthe 

35 

in the analysis of nonlinear wave interaction in a cold plasma. 

The averaged Lagrangian is first obtained in Section 4.1 for waves pro- 
pagating at an oblique angle to the static magnetic field. It was 
demonstrated in Chapter III that action is conserved in the linear theory. 
When nonlinear interaction occurs, action is exchanged among waves due 
to the wave coupling. The action transfer equations and the coupled mode 
equations are derived in Section 4 . 2 . In Section 4 . 3 , these equations 
are specialized to the case of waves propagating nearly parallel to the 
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magnetic field. In Section 4.4, the interaction among three cyclotron 
harmonic waves is considered. Coupled mode equations are presented in 
explicit forms suitable for computation. 

4.1 Averaged Lagrangian Density 

Taking the slow variation of wave amplitudes due to the mode coupling 
into account, we may modify Equation (3.1) to the form 

= £ ~k e3C PC ±0 k (^ t )] ' 

a(r,t) = 2 A (er,et) exp[i0 (r,t)] , (4.1) 

where the slowly varying amplitudes, (er, v,et) and A^ (e^r, et), 
may be written in terms of the Fourier amplitudes of linear theory, 

JkW and as 

^(er^et) = <^( «*> «t) [£( v) + e £ k (v)] , 

A k (er,et) = 0^(«r^«t) a k . (4.2) 

Here and can be regarded as correction factors accounting for 

the nonlinear effects, and 

a -k = ’ &-k = £k ‘ 

Consideration of the Lagrangian density expanded to second order 
in the perturbations sufficed for the linear theory of Chapter III. In- 
formation on nonlinear plasma phenomena is contained in the terms higher 
than second order. For sufficiently weak plasma nonlinearity, the third 
order term in the Lagrangian is much larger than those of higher orders, 
so only the second and third order terms need be considered here. 

The second and the third order terms in the averaged Lagrangian 
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density are obtained by substituting Equation (4.1) into Equations (2.2 6 ) 
and ( 2 . 28 ), and taking the average over a scale much longer in time and 
space than the period of the fast varying quantities, but much smaller 
than the scale of the slow amplitude variation due to nonlinear interac- 
tion. All terms in vanish, except those satisfying the resonance 

condition 


e k + V + V = 0 ' 

l 

which implies the synchronism conditions 


(4.4) 


k+k / +k // = 0, 0 ^+ 0 ^/ + (jq^// = 


0 


(^.5) 


( 2) 

For notational convenience, we will define A, such that 

7 k 


X,(2) . a (2) + a (2> .(2) . a (2)* 

\ “ A k A -k 7 A -k " A k 


(4.6) 


Although is always a real quantity, is not necessarily real. To 
first order in e, the part of the averaged Lagrangian density governing 
a weakly nonlinear plasma, i.e., the sum of the second and third order 
terms, is given by 


*-EhJ 2 4 2) ♦•2X [£ 


,( 2 ) 


+ e 


k+k '+k"=0 


f 

°kVV A kkV e x P[i( e k + e k 7 +v)] , 


(4.7) 


where the A* s are given by 
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4 M Af 0 < (2) + a * (2) ’ 

A k (2) = I K,oSj 2 + I <%(&> ' 4 x £ - «4 • (Sk + x x \) < 


« F <2) _ i , | e |2 . J_ |„ |2 

A k " 2 0 1 k> 2[1q k • * 


Vk'k" 2 ^ D k,C^k^ ’ ik' X £k" + 1 2 4Jk /: ~ C?k + ~ X 


In Equation (4.7)> the term involving j3 k has disappeared due to the 
first order equation of motion. 

4.2 Action Transfer and Coupled Mode Equations 

Variations of Z in Equation (4.7) with respect to 0 . 0 , and 0 „ 

K JK K 

give 


ft (kl\) + fr 

= It (mV) + 


•(^Kl\) - ft (MV) + 

h ■ (v'K*lV) - ‘(vw 



(M) 


where the summation is understood to contain all permutations of k^k'^k". 

O 

The first two equations are the so-called Manley-Rowe relations, 
giving information about the relative rates of transfer of action den- 
sity "N^j between the waves. The third equation gives the direc- 

tion and the absolute value of the rate of transfer. 

Variation of £ with respect to a. gives 

— K 
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(4.10) 


1N k(lt + z* 


I?) s - a -k'“-k" E 

k,k',k" 


A -k-k'- k " 


Similar equations for the k / and k^-waves are derived from the ct , / 

-k 

and a_ k // variations of X, respectively. Equation (4.10) is the 

coupled mode equation. It should be noted that the action transfer 

equation [Equation (4.9)] can be derived from the coupled mode equation 

by using Equation (3.12). 

( 3] 

Since N k , v^ and A_ k _ k /_ k /' are known for any wave, Equation (4.10) 

f ■*) 

can be solved in principle. The expressions for N. , v and A\ ' / ,/ 

k' k -k, -k ,-k 

in terms of field quantities are extremely difficult to evaluate in 
general, as indicated in Chapter III. A simpler form of Equation (4.10) 
is obtained in the limiting cases of waves propagating either nearly 
parallel or exactly perpendicular to the static magnetic field. These 
will be considered in the following sections. 


4.3 Quasiparallel Propagation (k « k ) 

j. z 

Substitution of Equations ( 3 . 36 ), ( 3 . 45 ) and ( 4 . 8 ), and some 
tedious but straightforward algebra, reduces Equation (4.10) to the 
following form first order in k^, 


IN, 


(it + ^k • h) 




= a -k /a -k* 


Z* Z* Z* R L* R* Z* L R* L* Z* 

//6 k VV +Q kk'k" e k e k'V + Q kk'k" e k e k 'V 


E (4v 

k,k / ,k // 

RR L* R* R* LL R* L* L* RZ L* Z* Z* LZ R* Z* Z*\ 

Q kkV e k VV + 9 kk'k" e k e k'V + s kk'k" e k e k'V + «kk'k" e k VV) 


(k.ll) 


where the Q’ s are given by 
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Q kkV 
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Q -k-k'-k" 


i. o a k Ji- 
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./ 2 
k v 
z x 
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2 X-i vVo 2a V 
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e O q .. 2 

2iT k z CU N \ 2 2 

^oA/O, 
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+ — k k v 
2 z z x 


1 + J-W 


%'o V;-i y 


i i 
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A/O^oA"-! W- 1 %0%1%'oVVl 


— 1 (— — -) 

v^,i^/oV;-i \v;-i “k,i / 
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+ — k k v 


2 z z x 


1 

i 1 + 1 \ + 1 

f 1 + 1 . \ 

A/oAaV;-! 

^'-1 “ k "°/ “mVxA 

Kl %o/ 
1 . 


1 + 1 


v\,A;oV;-i Wo %i / W- 1 a>k " 0 /J 


(4.12) 
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For illustrative purposes, the synchronism relations [Equation (4.5)] 
may be written as 


k=k+k. , cd = cd + cd. . 
~p ~s ~i ’ p s i 


(^• 13 ) 


In conventional parametric amplification terminology, the highest fre- 
quency component is called the 'pump' wave, and the lower frequencies 
the 'signal' and 'idler' waves. ^ With the use of Equation (4.13), the 
coupled mode equations for all possible cases can then be written in the 
form 


I 


U 

P 



s 


I 


W 

i 



+ 


U 

v 

~p 


+ 


V 


v 


w 


+ V. 



„u * 
e p = c uvw 

e s = C uvw 

e w - c 
l “ uvw 


V w 
£ £ , 
s 1 ’ 






(4.14) 


from Equation (4.11), where U, V and W can be either R, L or P; the 
coupling coefficient C T __„ is the proper linear combination of Q' s 
in the RHS of Equation (4.1l); v^ (U = R,L,P) is the group velocity of 
the U-wave, and and the slowly varying electric field, (!?, are 

K K 

defined by 



„U 


( er, et) 


a k ( er, et) 



(U = R,L or P; 
k = s,i or p) . 

(4.15) 


Here e must be replaced by e when it represents the electron 

K K 

plasma wave. The explicit expression for for right-hand polarized, 

left-hand polarized, and plasma waves can be obtained from Equations 
( 3 . 28 ) and (3*44) to yield 
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(4.16) 

Equation (4.1l) indicates that nonlinear interaction is possible (i) 
among three plasma waves, (ii) among one plasma wave and two circularly- 
polarized waves, (iii) among three circularly-polarized waves, and 
(iv) among two plasma waves and one circularly-polarized wave. The 
problem now remaining is to find the coupling coefficient in the 
coupled mode equation [Equation (4.l4)] for each of these cases: 

Case (i) (C ppp ): To first order in k^, all terms but the first 

one on the RHS of Equation (4.11) can be dropped in this case. The 
coupling coefficient then becomes 

C PPP = yi Q -psi > (^^T) 

-p,s,i 

z z 

where Q . is Q / „ in Equation (4.12), with k = -p, k' = s, and 

“P&l KK K 

, // . z 

k =i. Since Q^k 7 k 7/ c * oes not involve k^, the coupling coefficient 
remains unchanged for waves propagating exactly parallel to the mag- 
netic field. 

Case (n) (C^^p, C RPR ^ C LLP ; C LFL* ^RL* C PLR^ : In thlS case ’ only 
the second and third terms on the RHS of Equation (4.1l) contribute to 
the coupling coefficient, which becomes 


I R = I 


L* 




2 

%-l 


00 + 

C 


(Vx> 


2 -,\ 1 


r k - 2 u o 4 


n,o- 


W 15 ' 5 ' 1 ) “ C T.T.p( -p > -S > -1 ) - ( Q _r. = i “ ) } (^-iS) 


RRP 


LLP 


psi — sp — 1 


C PRL ^s-p + Q si-p^ ' 


(4.3-9) 
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'RPR 


C LpL and C pLR are obtained from C RRp , and C pRL , respectively, 


by interchanging s and i. As in Case (i), the coupling coefficients 


are the same for k 


0 and k = 0. 
l 


Case ( iii) ( C RRR > C LLL > C R RI/ ^R* C LLR ; C LRL^ : 


C is not of 
LLL 


interest because the topology of the dispersion characteristics does not 


allow the synchronism conditions of Equation (4.5) to be satisfied. 

z 

To obtain the coupling coefficients, the fi- * s have to be eliminated 
from the RHS of Equation (4.1l) by use of the small signal wave equation 
[Equations (3-40) and (3* 47)]/ 


,Z 1 , ZR R ZL L. 

’k = “ “ZZ ®k + A k V • 


(4.20) 


When this is done, the coupling coefficients C RRR , C RRL and C^ R are 
found to be 


RRR 


s.i 


RZ 

Q RR - _1_ (q r + q L ) 
V -psi ZZ - p s i “s-pl' 

i 


(4.21) 


C RRL ( ^ Sji) = C LR (_P>_S ^ i) = \ ^ 


-P;i 


■ 2 R, 

Q LL - - 2 - (q l + q r ) 

w si-p .ZZ v si-p ^is-p' 


(4.22) 


where A^ Z , A R Z and A ZR are given by Equation (3-47). The coupling 
coefficients C and C LRL are obtained from C RRL and C LLR , respec- 
tively, by interchanging s and i. The coupling vanishes as k 0 
in this case. 

Case fiv] (C , C , C , C , C„.,, C„„): All terms on the 
idse yxv; l^ R pp? PRP ; PPR* FliP* PPL 7 

RHS of Equation (4.11) must be changed into the form involving two 
£ Z, s with the use of the wave equation [Equation (3.40)]. This yields 
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RPP V ’ ’ ’ 


C LPP ( - p >- S ’- i) 


IE 


S,1 


a ZR 
RZ 

-psi ' RR v "'-psi 's-pi 
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(q r . + q l .) - 

v -nsi S-Dl ' 


ZR 

!e_ Q z 

ZZ ^si-p 


, (4.23) 




* . .V 

C ( -p, -s, - 1 ) 
PLP V ’ ’ ' 


E 


-p>i 


LZ i lr u ti \ _2_ 

Q si-p “ “LL (Q si-p Q is-p' A ZZ 


a zl 


A. 

1 


RZ 

s Z 
— 0 
ZZ i-ps 


, (4.24) 


where A RR , A LL , A RZ , A ZR and A ZL are given by Equation (3.47). We 
k ^ k k k K 

may obtain C ppR and C ppL from C pRp and C^, raapectively, by inter- 
changing s and i. The coupling vanishes as -* 0 . 


4.4 Perpendicular Propagation (k = 0) 

z 

As illustrated in Chapter III, for propagation oblique to the static 
magnetic field the polarizations of normal modes are extremely compli- 
cated. We shall study the relatively simple case of an electron velo- 
city distribution even in v . In this situation, ej^ and e^ can be 
considered as normal modes. These were introduced and termed the longi- 
tudinal and ordinary cyclotron harmonic waves, respectively, in Section 
3.5.2. As remarked at the end of that subsection, thermal effects on 
the propagation of the ordinary cyclotron harmonic wave are generally 
negligible. Therefore, for the coupling coefficients for the interaction 
involving the ordinary cyclotron harmonic wave, we shall approximate 
the wave by the ordinary mode for a cold plasma (see end of Section 4.5) • 
Since the phase velocity of the ordinary wave is much greater than 
the particle velocity, Equation (3.48) becomes 
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(4.25) 


-II _ G llll e ll ^ _ o-LII e " F Z - G ZZ e Z 

h. ~ G k \ ’ 5 k - G k e k » \ - G k 6 k > 


where the G* s are given by 


r ii ii _ a P n g-mi _ i a P n G zz _ _sl 

G k - m P -l,l » °k “ 1 m P 0,l,-1 * G k 


(4.26) 




Equation (4.10) is then written as 


/d 9 \ \ s / ii u* ||* ||* .0 ||* 0* 0*\ 

iN k at^'^K^-k^-k" 2 ^ KkV e k VV + Q kk'k" e k w) ’ 

' t.tt ' 


k,k,'k 


(4.27) 


where the superscript Z has changed to 0 to represent the ordinary 
cyclotron harmonic wave, and the Q' s are given by 

iv - - 1 v k l ([ G “ “ s (Mk">-< ““(W] 

X jo”, cos^,-^,,) - G k [, sln( 4> k )J\ , 

iri k l -ir— ([ D -k,o +i K'-v)][ G " 

m n n 1 - 


“k'v 


(4.28) 


Substituting Equations (4.5) and (4.26), using the Bessel function iden- 
tities of Equation (3.42), and noting that 


( D k, 0 _iai k,n) ex PU[»(H s ) + X k Sin ^ _< ^k^^ = 0 * 



(<M k ) + X k /Sin(<H> k ) ] }d<() = 2nJ n ( ) exp( in<j> k „ ) , 


(4.29) 
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we can reduce Equation (4.28) to the form 




. “O' 1 . // 2 , . v 


_ ir » < J »(V) J n^k-) J M .n(V)> 

^^m^/n 


' ex P[ i( ) ] ex P[ -i(W ) ] 


I exp[ - i( (j) k ,-(j) k // ) ] exp[ - i( <j> k ,-(jy ) ] ' 


^k/n+l 


^/n-l 


exp{i[m( <i) k -(J> k // )+n( <j> k ) ] } , 


.0 • 6 0 q v 

Q kkV = 1 w k 


j_ v Ln 4l\ 

<\ 7 \,n 


/ ex P[ ) ] exp ( -i ^k -< l ) k // ^ ^ ' 

X + 


^kjn+l ^^n-l 


(4.30) 


As in Equation (4.10), the summation in Equation (4.27) contains all 
permutations of k, k', k /7 , while those in Equation (4.30) are double 
summations of m and n over all positive and negative integers. 

Equation (4.27) can be written in the form of Equation (4.l4) with 



(4.31) 


These result form substitution of Equations (3.28) and (3.44) 
into Equation (4.15). Possible interactions are, from Equation (4.27)? 
(i) among three longitudinal cyclotron harmonic waves, (ii) among one 
longitudinal and two ordinary cyclotron harmonic waves. The coupling 
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coefficient for each of these cases can readily be obtained from Equa- 
tion (4.27): 

Case (i) (C ): In this case, only the first term on the RHS of 

Equation (4.27) contributes to the coupling coefficient, which becomes 



= X) q -p is ’ 


(4.32) 


where Q 11 . is qJJ, / // with k = -p, k' = s, k" = i. 

“ p 1 S KK K 

Case (ii) (C C C ^ ): In this case, only the second term 

' ' v ||OCr 0)| 0 ^ 00 || ' 

on the RHS of Equation (4.27) contributes to the coupling coefficient, 


which becomes 



(4.33) 


c = N ' Q°. (4 

0(1 0 / J *81 -p 

-P>i 

We can obtain C from C by interchanging s and i. 

UUJI U|| u 

The infinite series forms of the Q' s in Equation (4.28) are 
impractical for computation. Alternative expressions for them are 
obtained in Appendix B. 

4.5 Discussion 

We have been following essentially the averaged Lagrangian method 
developed in cold plasma wave theory. Since all variables are de- 
fined in phase space, the method is modified for warm plasma by 
adding the correction factor j^, which accounts for nonlinear 
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effects on velocity. However, it has not produced any additional 
Euler-Lagrange equations because the term involving j3^ has dis- 
appeared from the averaged Lagrangian in Equation (4.7). 

General solutions of the coupled mode equations [Equation (4.l4)] 

4 

are complicated, but can be obtained in terms of elliptic integrals. 

For the special case where one wave is considerably stronger than the 

other two waves, however, solutions can be obtained in a simple form. 

If the pump wave is very much stronger than the signal and idler, then 

U V W U 

the rate of change of £ is of order & £. , and £ may be assumed 

P si p 

to be constant . The signal and idler waves are determined by the 
equation 


[(It ♦£•§;)( 


a . w 


d 

dr 


) - = 0 - < lt -35) 


/ X A . w f|_ „U,2 /T V T W*a/2 

where the growth (or decay) rate, y, is given by ^l C uvvrp' ' s*i 

V W . . 

The electric fields £ g and £^, are either growing or decaying, depending 

on the sign of I g I^ . From Equations (3.22) and (4.15)? this sign is 

V W V W 

the same as that of [a) co. J( K. 1 where H and H. are the energies of 

L S 1 S 1 J S 1 

V W 

the signal and idler waves, i.e., £ g and £ ± are exponentially growing 
(decaying) when [oyixMgK^] >0 (< 0) . The growth (or decay) rate can 
be calculated from Equation (4.35) for a specific unperturbed electron 
velocity distribution, and wave numbers and frequencies satisfying the 
synchronism conditions. 

Table 4.1 shows all possible nonlinear interactions among normal 
modes propagating nearly parallel or precisely perpendicular to the 
static magnetic field. The third column of Table (a) indicates the 
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Table 4.1 

* 

All Possible Nonlinear Interactions Among Normal Modes 


(a) Quasiparallel Propagation 


Involving Waves 

Pump-Signal 
and Idler 

Coupling 

for k -=0 
l 

Previous Investigation 

— 

microscopic 

treatment 

cold plasma or 

hydrodynamic 

approximation 

(i) Three plasma waves 

P P,P 

yes 

Refs.28,33,34 

Refs. 17,20 

(ii) One plasma wave 

and two circularly— 
polarized waves 

R R, P 

L L,P 

P R,L 

yes 

Refs. 33,34 

Ref. 26 

(iii) Three circularly- 
polarized waves 

R R,R 

no 

Ref. 34 

Ref. 27 

l , 

no 

None 

None 

(iv) Two plasma waves 

and one circularly- 
polarized wave 

R P,P 

L P,P 

P R, P 

P L, P 


( b) Perpendicular Propagation 


Involving Waves 

Pump-Signal 
and Idlerf 

Previous Investigation 

microscopic 

treatment 

cold plasma or 
hydrodynamic 
approximat ion 

(i) Three longitudinal 
cyclotron harmonic 
waves 

II II > II 

Refs. 32,80 

Ref. 25 

(ii) One longitudinal 
and two ordinary 
cyclotron harmonic 
waves 

0 || ,0 

Refs. 30,65 

Ref. 22 

■n 

None 

Ref. 22 


(For the electron velocity distribution function even in v z . Inter- 
actions between parallel propagating waves and perpendicularly pro- 
pagating waves are not included.) 

+ 

( || : Longitudinal cyclotron harmonic wave - extraordinary wave in 
cold plasmas, 

0: Ordinary cyclotron harmonic wave - ordinary wave in cold plasmas.) 


68 








































interaction in the limit of k =0. As shown in Section 4.3* when 
three plasma waves or one plasma wave and two circularly-polarized 
waves are involved, interactions still occur for all waves propagating 
exactly parallel to the static magnetic field, and their coupling co- 
efficients are the same as those for quasiparallel propagation. The 
remainder of the interactions in the table vanish for exact parallel 


propagation because their coupling coefficients are first order in . 

The last two columns indicate the interactions previously inves- 
tigated by other methods for cold or warm plasmas. The coupling co- 
efficients obtained in Section 4.3 for the first five situations in 

34 

Table (a) are identical to those derived by Kim, Harker and Crawford. 

The first case in Table (b) involving three cyclotron harmonic waves 

32 80 

was first investigated by Harker and Crawford" 5 , and later by Wagers. 

The former considered collinearly propagating waves (<j) = <(>/ = <(>,//). 

The latter generalized the analysis to allow for waves propagating at 
arbitrary angles in the plane normal to the magnetic field ( (j) ^ (|) /*£ (j) ,,) . 

K K K 

The coupling coefficient C in Equation (4.32) is in a form different 

II II II 

from those obtained by these authors. While they could not derive the 

Manley-Rowe relations analytically from the coupled mode equations, the 

present coefficient is written in a form symmetric in -p,s,i and 

automatically satisfies the relationship; a useful indication of the 

power of the Lagrangian method. The coupling coefficient C _ has 

°!l° 

been derived by Boyd and Turner ^ for a Maxwellian velocity distribution. 
Their result follows readily from Equation (4.34). 

Table 4.1 also shows the new interactions obtained in this chapter. 
Those involving three circularly-polarized waves are the combinations 
of R-R, L and L-R, L. Typical synchronism conditions for these inter- 


actions are shown for a cold plasma in Figure 4.1. For a warm plasma, 
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thermal effects make only a very slight difference to the shape of 
the upper branch of the R-wave and the L-wave branch. The lower 
branch of the R-wave can be strongly modified, and the three inter- 
acting modes in Figure 4.1(a) can all be in the first quadrant of the 

co- k plane, as shown in Figure 4.1(b) 

As shown in Table 4.1, the interaction involving two plasma waves 
and one circularly-polarized wave has four cases: R-P,P; L-P,P; P-R,P; 
and P-L,P. The topology of synchronism conditions of the L-P,P case is 
illustrated for a cold plasma in Figure 4 .2( a) . Figures 4.2(b) and (c) 
shows synchronism conditions of the R-P,P case for a cold plasma. When 
oiy < 20^ (or < 30^/2), the right-hand polarized wave is on the 
upper branch while it is on the lower branch when cn c > 20^. As men- 
tioned above, the lower branch of the R-wave can be strongly modified 
by thermal effects and the interacting R-wave can always be on the 
lower branch. The last two cases (P-R,P and P-L,P) have no solution 
in cold plasma theory since the synchronism conditions obviously cannot 
be satisfied. For a warm plasma, the whistler and the plasma waves 
can be unstable for a nonMaxwellian-velocity distribution. The possi- 
bility also exists of these waves propagating with either positive 
or negative small-signal energy. It would therefore be of great 
interest to make further detailed studies of the R-P,P and P-R,P 
cases to determine the relative growth rates of the small-signal and 
nonlinear instabilities, and to see which should dominate in experi- 
mentally relizable situations. 

A new interaction of perpendicularly propagating waves treated 
here is the case ||-0,0. Figure 4.3(a) shows its synchronism condi- 
tions for the 0-wave representing the ordinary wave in a cold plasma. 
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Figure k.2 Synchronism conditions for interaction among one 
circularly-polarized wave and two plasma waves. 

(a) Pump L-signal and idler P,P; (b) Pump R-signal 
and idler P, P for o^. < 3%/2; (c) Pump R-signal 
and idler P,P for a> c > 2o^. 
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cold plasma; (b) O is the ordinary cyclo- 



Since the dispersion curve of the ordinary wave is inside the light 
cone, so is the pump wave as shown in the figure, where the ||-wave 
is non-dispersive . Unless the magnetic field and plasma were homo- 
geneous, synchronism would not be maintained. Furthermore, because 
of the zero group velocity, collisional effects exert a strong 
influence and may invalidate the theory. Therefore, this region of 
the y-wave is not of great interest. Figure 4.3(b) shows the 
synchronism conditions for the O-wave representing the ordinary cyclo- 
tron harmonic wave. In the region satisfying the synchronism con- 
ditions, the waves are dispersive. How good our approximation is for 
C [Equation (4.33)] open to doubt since we assumed cold plasma when 
simplifying it. The whole question is probably somewhat hypotehtical 
in practice; the interaction would be extremely difficult to detect 
experimentally because the propagation bands of the O-wave are very 

79 

narrow. ^ 
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Chapter V 


HYDRODYNAMIC APPROXIMATION 

The hydrodynamic approximation for plasma wave propagation is 
appropriate when the thermal velocities of the charged particles are not 
comparable to the wave phase velocity. The plasma can then be described 
in terms of macroscopic quantities such as its density, drift velocity, 
pressure, etc. These quantities are related to each other by moment 
equations obtained by carrying out the velocity integration of the Vlasov 
equation. If the adiabatic assumption is made, the infinite set of 
moment equations so obtained can be truncated at the third. The first, 
second and third moment equations, known as the continuity, momentum 
and state equations, respectively, together with Maxwell's equations, 
then provide a closed set of equations. 

It should be possible to describe the macroscopic properties of 
plasmas by Hamilton's principle by formulating a suitable hydrodynamic 
Lagrangian. In this chapter, such a Lagrangian will be derived by two 
different methods. In the first, the hydrodynamic Lagrangian density 
is obtained for the momentum equation and Maxwell's equations by inte- 
grating the exact microscopic Lagrangian of Chapter II in velocity 
space. In applying Hamilton's principle, the continuity and energy 
equations are then used as subsidiary conditions. To describe linear 
propagation and wave-wave interactions, this Lagrangian is expanded 
about the equilibrium state. Because of some mathematical difficulties, 
the expansion will be carried out here for only three illustrative 
cases: one-dimensional compression parallel to the static magnetic field, 
two-dimensional compression perpendicular to it, and three-dimensional 
compression. In the second method, the expanded hydrodynamic Lagrangian 


75 



is obtained directly from the expanded microscopic Lagrangian of 
Section 2.2. In contrast to the first method, this does not require 
any subsidiary condition, and avoids the mathematical difficulties 
involved in the expansion of the Lagrangian. 

In Sections 5.1 and 5.2, the Lagrangian is formulated by the 
first method. The second and third order Lagrangians are then obtained 
for the three illustrative cases mentioned above. The second method 
is presented in Section 5 . 3 . Section 5.4 describes how the ion 
contribution modifies the first order Euler-Lagrange equations and the 
coupled mode equations. In Section 5 . 6 , the nonlinear interaction of 
two electron plasma waves and an ion acoustic wave, and of an electron 
plasma wave, an ordinary wave and an ion acoustic wave , are studied as 
an application of the Lagrangian by the first method. The application 
of the second method is carried out in Section ^. 6 , where dispersion 
relations are obtained from the second order Lagrangian of the second 
method and are compared to those in Chapter III. As an example of 
nonlinear wave interaction, the interaction of an Alfven wave with two 
electron plasma waves is examined. 


5 .1 Formulation of the Lagrangian Density - Method I 

The hydrodynamic equations obtained by integrating the Vlasov equa- 
tion [Equation (2.6)] in velocity space with the adiabatic approxima- 

. . 81 

tion are 




u = 0 , 


(5-1) 


du a 
dt dr 




q 

m 


(E + u X B) , 


(5-2) 
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ds 

dt 


a ♦-> «-* 

— s x B + s 

ITl /v» ~ 



- transpose , 


(5-3) 


where the local plasma density, n, the pressure tensor, P, and the total 
time derivatives, d/dt, are respectively defined by 


n 



} 



d _ 3 d 

dt “ St ~ ' Sr ’ 


(5-4) 


and the tensor s, and the random velocity, w, are related to the pressure 
tensor, P, and the drift velocity, u, respectively by 



mn 


y 


w = v - u . 


(5-5) 


The diagonal elements of s are the squares of the thermal velocities 
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in each coordinate direction. Equations (5*1) “ (5-3) are the equa- 
tions of continuity, momentum, and state, respectively. Together with 
Maxwell's equations, they constitute a complete set of equations describing 
adiabatic processes in plasmas in terms of macroscopic quantities. A 
Lagrangian for the hydrodynamic approximation should include these equations. 

Analogously to integrating the Vlasov equation in velocity space in 
order to obtain the moment equations, we may attempt to integrate the 
microscopic Lagrangian in velocity space to obtain the Lagrangian for the 
hydrodynamic approximation. Transforming the Lagrangian coordinates 
( r 7 ,v ') to Eulerian coordinates (r,v), and integrating Equation (2.2) 
over the velocity, yields 


„ nm 2 , nm m <-* , 0 

Z = — u + — Trs + nq(cp - u • A) + — 

c_ c_ ^ cL 



(5.6) 


However, it should be remembered that the conservation law in phase 
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space was used as a subsidiary condition in connection with Equation 


(2.2). The information contained in the Vlasov equation has therefore 


been lost from this integrated Lagrangian. It has been shown in other 
Stanford work by Peng that changing the sign of the term Tr s is 
necessary to obtain the required Lagrangian, 

£ = u 2 - Tr 7 + nq(cp - u • a) + ~ ( Vc P + ^ ( 7 X A) 2 . (5-7) 


It is easily seen through the use of Equations (5*1) and (5*3) as subsi- 
diary conditions, that the variation of this form of £ with respect 
to r gives Equation ( 5 . 2 ). As in Chapter II, the variations of cp 
and A give the Maxwell equations . 


5 .2 Expansion of the Lagrangian Density 

To follow the procedure of Chapters III and IV for small-signal wave 
propagation and wave-wave interactions, it is necessary to expand the 
Lagrangian. As in Chapter II, two corresponding cells satisfy 


£ = £o + -l> 


n(r,t)d 3 r = n^Jd 3 ^ , 


(5.8) 


where n (r ) is the equilibrium density. For a small perturbation, 

\J 

the local drift velocity and the tensor s of the perturbed cell may 
be expanded about the cell in equilibrium (see Section 2.2) to give 


= U.C^q) + > I = AT + (£ • v n)i , 




5t 


0 ' 


(5-9) 


B^) = S,(*b) + + + £3^) + 


(5.10) 


where U is the drift velocity of the cell at and the subscripts 
on s indicate the order in §. The second relation in Equation (5*9) 
has been obtained by an argument similar to that given in Appendix A. 
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The expanded Lagrangian can be obtained from Equation (5*7) by ex “ 
panding Equation (2.l4) about in a Taylor series, and substituting 

it together with Equations (5-8) - (5-10). However, it is still re- 
quired that ^ be expressed in terms of ij, and a. In principle, 

this could be done by expanding the ten simultaneous differential 
equations expressed by Equations (5*1) “ (5*3) (one for n, three for 
u, and six for s) to the desired order in perturbation. In practice, 

M 

obtaining the solution for is extremely complicated. For s^, for 

example, the number of differential equations to be solved simultaneously 
after expansion is 26 (2 for n to the second order, 6 for u to the 
second order, and 18 for s to the third order) . The solutions for s 
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are very much simplified for the special cases considered by Delcroix, 
i.e., linear adiabatic compression parallel to the magnetic field, cylin- 
drical adiabatic compression perpendicular to the magnetic field, and 
spherical adiabatic compression. In what follows in this section, we 
shall restrict the discussion to these three illustrative cases. 

5.2.1 Linear Adiabatic Compression Parallel to the Magnetic Field 
Let the static magnetic field be uniform and directed along 
the z-axis; the velocity distribution at every point have rotational 
symmetry about the z-axis, and all perturbations be independent of 
x and y. Thus, 



s =s = s , s = s , s = s =s =0 
xx yy x ’ zz || xy yz zx 


(5.H) 


where s and s can be related to the electron thermal velocity, 

x || 
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v^, and temperature, T, by 


< 2 - 
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s 

1 


2 

kT 

ti ± 

2 “ 2m 


1 



(5-12) 


This case is the same as that of electrostatic waves in a uniform plasma 

♦4 

propagating parallel to the magnetic field. The six equations for s 
in Equation (5.3) reduce to two, 


ds 


dt 


= 0 , 


' n ' 


(5.13) 


Expansion of Equation ( 5 . 13 ) is carried out in Appendix C. 

Substitution of Equations ( 5 . 8 ) - (5.10), and (C. 8 ), into Equation 
(5*7) gives the second and third order Lagrangians as 


» _ *2 3 

£ 2 ~ 2 5 z 2 0 11 


i*h\ 2 
\* r / " V 


^ e + e 2 , 
^z z 2 z ’ 


£ = 2 n m s 

3 0 11 


/ay 3 _ v .2 

\dz / 2 Bz ’ 


(5.14) 


where e z is given by Equation ( 3 . 16 ), and the subscript 0 has been 

dropped from s . By summing £ and £ , the averaged Lagrangian can 

II ^ 

be written as Equation (4.7) with and A^V' ®iven by 


.(2) 0 2 | _z 1 2 3 I, _zi2 F z z* S 0 1 zi2 

A k " 2 %0 lil 2 n 0 mS n l k z ? kl “ n 0 Q ^k e k + 2~ l e J » 


- i 2 n o ms „ vx + 1 v k Ai<- ■ (5.15) 


where has been written for - k^ + i^cd, . Here it should be 

noted that $ = 0 implies U = 0, from the zeroth order of Equation 
( 5 . 2 ) for a uniform plasma in equilibrium. 
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5 . 2.2 Cylindrical Adiabatic Compression Perpendicular to the 
Magnetic Field 

For compression purely perpendicular to the magnetic field, 


we may put 


\ 


§ 


z 




s 

xz 


s = 0 

yz 


) 


B = B + b , B = B = 0 . ( 5 -l 6 ) 

z 0 z ’ x y i 

If we further assume a symmetric compression about the z-axis, we may 
also put 


s = s = 

XX yy 

s , s = 

L ZZ 

= s , 
II 

s = 0 . 

xy 

(5.17) 

The six equations for s 

in Equation 

(5-3) 

then reduce to 



s \ 

ds 



Sf( 

0 

11 

II 

dt 

11 

0 

s* 

(5.18) 

He 

dx 

35 

" 3y ' 

35 

^x 

w 

3? 

JL 

9 x 

(5.19) 


Expansion of Equation (5.18) is carried out in Appendix C. 

Of Equations (5.16) - ( 5 . 19 )> we only need Equation (5.18) for the 
expansion of the Lagrangian in Equation (5.7). Let us consider the 
implications of Equations (5.16) and (5.I9) . In Equation (5.16), the 
second relation implies that the waves under consideration are propa- 
gating perpendicular to the magnetic field. Furthermore, since the 
time-varying magnetic field of the wave is in the z-direction, it can 
be seen from Equation (3.16) that e^ is polarized in the xy-plane. 

It follows that there is no time varying displacement in the z-direction. 
Equation (5. 19) expresses the rate of strain relations,^ which become 
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3^/3x = SI- /3y an< 3 3^ x /3y = - 3§y/3x for the steady state. They 

imply that the motion due to a wave is of uniform dilatation in 

86 

the xy-plane with rotation about the z-axis. This motion is only 
possible for a wave whose frequency is well below the cyclotron fre- 
quency, so that the displacement of a fluid element due to the wave 
is sufficiently randomized by the gyration for symmetrical compression 
to be approached. A wave satisfying both Equations (5.16) and (5.19) 
is the compressional Alfv4n wave with a frequency well below the ion 
cyclotron frequency, and propagating perpendicular to the static magne- 
tic field. 

Although compressional Alfven waves propagating perpendicular to 
the magnetic field can satisfy the synchronism conditions for wave-wave 
interaction among themselves, we shall only consider the second order 
Lagrangian here. Substituting Equations (5-8) - (5-10), and (C.10), into 
Equation (5.7) gives the second order Lagrangian as 


Y . 

2 2 


F 2 - n ms ( V • F ) 2 - 

b X 0 x [} X ^x' 


0 3§ 3? 3§ 35 

2 ’x ^y , ^x 3 


3x dy 


_Y| 

3y 3x J 


V ' S 0 2 1 2 

n q F . e + -X- F • F x B + -^ e - ~ b^ , 

0 --3i — 1 2 '-Sx ~ 2 2^^ z * 


(5-20) 


where the subscript 0 has been dropped from B and s . The averaged 
Lagrangian then becomes 



n„ms 
0 x 




2 





(5.21) 
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5.2.3 Spherical Adiabatic Compression 


In this case, even though perturbations occur, the plasma 

4-4 

preserves isotropic properties. The elements of s can therefore be 
written as 


s = s = s = s . s = s = s =0. 
xx yy zz ’ xy yz zx 


( 5 - 22 ) 


Equation (5-3) then yields 



dx dy 



> 



(5.23) 

(5-24) 


If the first relation of Equation (5.24) is identically zero, then 
Equation (5.24) describes a combined translation and rotation as in 

the case of rigid bodies. The first relation is not zero in general. 

86 

This implies the addition of a uniform extension (or contraction) 

Q r 7 

to the translational and rotational motion. Fluid motion of a uni- 
form extension (or contraction) is possible in a plasma in which 
elastic collisions conserve energy and momentum, and in which the 
frequency of the wave is much lower than that of collisions. The 
motion of the fluid elements due to the wave is then randomized by 
collisions during a period of the wave. The translational and rota- 
tional motion is due to the static electric and magnetic fields. 

Substitution of Equations ( 5 . 8 ) - (5.10), and (C.ll), into Equa- 
tion (5.7) gives the second and third order Lagrangians as 
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JU = 


n m , „ n q . 

-§-^-?V<’-S 2 + -§-S 


rs» 


§ X B 


n 0 q l ■ <£ + £ x i) * t e2 ' kr 0 b2 • 


*3 - ' fl ”o ms (’ ■ !> 3 + 5 ■> 0 "=(15:W)( q ■ ?) 


n o q n o q • 

--H-^v(e+uxb)+-H-?x( 5 xb) , 

rv-/ <~w- rs-/ ^ rv rv»> 


(5-25) 


where the subscript 0 has again been dropped from s, and in obtaining 
the first two terms in from Equation (C.ll), the total divergence 

terms are neglected. The averaged Lagrangians are then written as 
Equation (4.7) with 


(2) V 2 


“k,0^ f V s 4 ' £l +1 2 “i,0 4; ■ -5k X £ 


' n 0 q 4 ' ^ik + ~ x -V + 2° l^k |2 ■ 2j^likl 2 > 

4kV' - - 1 f; V s + 1 1 V s Kfe ; &Jk'V 

+ 1- |-4'ik" ; i.(& + £ >< £i 1 ) - 1 -§-%olk x (4' ,< V>-(5-26) 


5 • 3 Formulation of the Lagrangian Density - Method II 

The hydrodynamic Lagrangian in Section 5 .1 is formulated from the 
exact microscopic Lagrangian expressed by Equation (2.2) . Because of 
the mathematical difficulties involved in the expansion, however, its 
practical use is restricted to very special cases. In order to avoid 
these difficulties, we will obtain the expanded Lagrangians for the 
hydrodynamic approximation directly from the expanded microscopic 
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Lagrangian of Equations (2.2 6 ) and (2.28). 

Let the random velocity at equilibrium be W. The velocity v- 

rv 

coordinate system in Equations (2.26) and (2.28) may be changed to the 
random velocity W-coordinate system via 

v = u + w , (5.27) 

rs-* yvy 

where the subscript 0 is dropped from the equilibrium velocity, v^. 

This change of coordinate system results in the following transformation: 


S_ 

3t St * 



s_ s_ 

Sv "* 3W • 


(5-28) 


For a plasma uniform at equilibrium, the second term on the RHS of the 
second transformation in Equation ( 5 . 28 ) can be dropped. In the absence 
of a static electric field, the time differential operator, D, defined 
in Equation (2.19)? may be written as 

D = Ht - CS*“c> ■ (Z- 7 > > 3t*ft + D zfe- (5-29) 

Here, the subscript 0 is dropped again from t, r and V. Let the plasma 
temperature be very low, i.e., the mean value of the Larmor radius, and 
the mean thermal velocity of particles are very small compared to the 
wavelength and phase velocity, respectively. Since the spread of the 
velocity distribution function is small for this situation, the displace- 
ment ^ may be expanded about v = U as 

rvy y^y 


5 (u + w) = 5 (u) + (w 


yv-y yvy 


/Ny y^y 


V IW + k W 1W + 


(5.30) 


Substituting Equations (5.27) - ( 5 . 30 ), dropping the terms of 7 V E 

w 

and those of higher order than second in W, and using /w>= 0 , 
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reduces Equations (2.26) and (2.28), to 


-"o<! 


d§ 2 d§ 


+ i* • •§* (%?> 7 *s 


* ■ 3? • (™ OT w)l - »(* • VS ■ * x i) + iO™ :7 VS' £ x 5. 


+ §(I' 7 S + £' 3t (''wD-.**®- V 2 -23T b2 ' <5- 31) 


n o q / d J 

= • 5 X b - §§:*(e + u X b) 


+ [W • 7 ? + w. — ( 7 „?) - — W X B • V „§) 1 • (W . 7 j?X b+ (WW: 77 J § • § Xb 
L ~ .2. ~ dt v W~ m ~ ~ W-2/ J ~ W~ ~ Wr.o ,2, ~ 


(I- V I<S- VS ;7 ^ + £ x t)-<»- 7 w>(SS> :7 (* >< £)) • (5.32) 


Though it is not a necessary restriction, the velocity distribution 
is now assumed symmetric for the sake of simplicity, so that the off- 
diagonal elements of s^ = vanish. The averaged Lagrangian 

then yields Equation (4.7) with 


A (2) _ x (2) + (2) + (2) (2) 

A k - k 0 A 00 A 10 X 11 ' 


( 5 - 33 ) 


A ( 3 ) _ ,( 3 ) + ,( 3 ) , ,( 3 ) , .( 3 ) , _ ( 

A kk 'k “ X 0 + X 00 + X 10 + X 11 + C -‘ 


( 5 - 34 ) 


The X V 1 s in Equations ( 5 * 33 ) are given by 
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x o 2) -tv“k,o I 4 I 2 - V4-<^ + S x ik ) + l n o q< \;c4- <^ x 5) 


+ i'oW 2 -4i4i 2 . 


x 0 2) = I n o‘ k a s j4l 2 ' 


^10 


,( 2 ) _ 1 
ni ~ 2 


iqn G^,o # cpY B c«Pp ; 

v 2 54 

2m e opY S QD! 0 9W 


* 

* 

’aa£k ‘ aw 
a 

- ln o“4 ' 1 

3 Jk / 3 J 

sT + ^! 

t aw ‘lav 

! 1 2m 

a. \ 

a/ 




e^ v s B„ 


?v. 

3 ?k 


P * . 
2 k ‘ SW 

a 


aw. 


^ ^ R 

)S on aw a • &*„ * ~ 


/ a 4\* „ « 

• ^3«r) * 2 . - n o q€ apv s oD£ « \3wj ' 


(5-35) 


where a,p, and Y represent in any order the variables x, y, and z, the 

subscript 0 has been dropped from s aa ’ and the property of the skew- 

symmetric tensor, e .... that a . bxc = e a b e has been used. 

Op Y ~ ~ ~ OCp y a p Y 

The rules for the subscripts and superscripts are as follows: the 

subscripts on k and zero older variables indicate the components of 
the coordinates; the subscripts on o> and first order variables indi- 
cate the participating wave vectors. The coordinate components of 
these variables are indicated by superscripts. The only exceptions 
are || and x, the use of which was explained in Section 3*5 
The X’s in Equation ( 5 . 3 ^) are given by 
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! 



x o 3) - ; "o q “t,olk • ik' x £k" + | n o*(^k + E x Ek > ! W > 


,(3) _ o 
A oo 0 > 


4? - - f n o qs c« (k a- k a ) h ' 5T" * V + f VE ' e opY s cn b k 5 k' 

a a 

i k , «<- 

* 2 “o q E • ik ' e QpY s aA 55- ' 

a 

(3) i 3 Jk S ik' n 0 q2 3 lk' 

X ii = 2 n O q %O s oa W ‘ aw~ x 5k" ~ ~&r e apY S co B p dw" ’ Bw~ x 5k" 

a a y a 


i 3 Jk' s Jk" 

+ 2 V s QQK'^~k + £ x 5k ^ : sw aw~ 

a a 


5.4 Inclusion of Ion Motion 

Since the second and third order Lagrangians have now been obtained 
in Sections 5.2 and 5-3? we are in a position to follow the procedure 
described in Chapters III and IV to obtain linear dispersion relations 
and the coupled mode equations. Before doing this, however, since ions 
play an important role for the hydromagnetic waves, we must consider 
how the inclusion of ion motion modifies the Lagrangians and the Euler- 
Lagrange equations. 

Including the ion motion, the Lagrangian may be written as 

£ = / e + Z 1 + / , (5.37) 

Pe 

where the electron contribution, £ > is written for the terms in the 

Lagrangians in the previous sections which contain the displacement 
vector §. In order to denote the electron contribution, superscripts 
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(e) must be appended to m, u, U, ai , a) , s, ^ and ^5 of Sections 

^ ^ t M* ^ ^ ^ 

Pi P 

5.2 and 5 . 3 . The ion contribution, jf , is obtained from £ by 
changing the superscripts (e) on these variables to ( i) with the following 
changes in sign: 


q -* -q , 




( e ) ( i) 

01 ^ 01 


( 5 . 38 ) 


In the linear theory, the equations for the ion motion to be derived 
from the variations in and V are therefore the same as 

those for the electron motion, except for the superscripts and sign 
changes just noted. In Maxwell’s equations, which are obtained by varia- 
tion of the fields, the current and charge densities consist of the 
contributions of both electrons and ions, i.e., 


( e) . ( i) ( e) ( i) 

= ' + ;T ' , P = P v ' + P v ’ , 




where j^^ and p^"^ are respectively obtained from j^ 6 ^ and p 


(5-39) 

(•) 


with the same changes as noted for S ?* . The action, N^, and Hamiltonian, 
K^, should also be written in analogous form to Equation (5-37). 

In the nonlinear theory, the correction factor Q^, defined in 
Chapter IV, is the same for both electron and ion motions. The other 
nonlinear correction factor, ( 3 , does not appear in the hydrodynamic 
approximation because £ and are no longer functions of v. The 

rZ/ Vr*-' 

Lagrangian, the action transfer equation, and the coupled mode equation 


can then be written as Equations (A. 7), (^*9) and ( -U . 10 ) with N k , 
(3) 

and A^/j^// including the effects of ion motion. 


( 2 ) 


5 .5 Application of the Lagrangian by Method I 

As an application of the Lagrangian formulated by the first method, 
we shall consider the interaction of waves which belong to the parallel 
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compression case treated in Section 5-2.1. Cylindrical compression 
need not be studied since the waves do not satisfy the synchronism 
conditions. The case of spherical compression is for collision- 
dominated plasmas, and is outside the scope of our work. It is 

Oo 

under extensive investigation at Stanford by Peng. J 

The waves in the parallel compression category are electrostatic 
waves: electron plasma and ion acoustic waves. In this subsection, we 
shall obtain the linear dispersion relations and the coupled mode equa- 
tions for these waves. The ordinary mode does not belong to this cate- 
gory since it has variations in the xy-plane (k^ 4 0 ) , and has a time 
varying magnetic field. It thus violates the assumptions of Equation 
(5.11). However, since the wavelength of the ordinary wave is usually 
much longer than those of electron plasma and ion acoustic waves, as 
shown in Figure 5 . 1 , the dipole approximation for this wave may be 
satisfactory when it interacts with waves propagating axially. Con- 
sequently, we shall also consider interaction involving the ordinary 
wave . 

5-5-1 Linear Theory 

For small-signal propagation, the terms of order 0(e) in 
Equation (4.7) can be neglected. Since only a single wave number, k, 
is necessary, the Lagrangian may be written as 





(5- 1 to) 


where A 
respect 



is given by Equation ( 5 . 15 ). Variations of 

_z(e) z(e) 

and a ^ ' give, respectively, 



with 
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kp«ks,kj 



(b) 


Figure 5.1 Synchronism parallelograms for interaction of one ordinary 
wave (pump), one electron plasma wave (signal), and one 
ion acoustic wave (idler). 
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(5.41) 


E z ( e ) 




2 

3 k s 
z 



( e ) 

4,0 H ^ 


k l/ 6 ) 
z z 


•T Z ( e ) J. (*) 

-lk e«e = p v + p v 
z 0 k r 


.(«) 


= -‘‘Vz’k 


F Z ( e ) 


The first expression in Equation (5.4l) is the equation of electron fluid 
motion. A similar relation for ions is obtained by the proper changes 
in Equation ( 5 . 38 ). The third expression is Poisson's equation. Com- 
bining the relations in Equation ( 5 . 1}-1 ) gives the dispersion relation as 


= 1 


-E 


(a) 2 

4 


a 


(a)2 

4,6 


- 3 k 2 s(“> 

Z II 


= 0 , (a = e,i) 


(5.42) 


where is the equivalent plasma permittivity. This is the familiar 

dispersion relation for an electrostatic wave. A simpler form of dis- 
persion relation for an electron plasma wave can be obtained by assuming 
the frequency and the phase velocity of the wave to be very high com- 
pared to the ion plasma frequency and the average ion thermal velocity, 
respectively, i.e., 



» cu, 


(i) 


N 


and 


z 11 


( 5 . 43 ) 


Equation (5-42) then reduces to 


( e ) 2 (e )2 . -.2 (e) 

“k ,0 * % + 3“z s ’ 


.(«) 


kT 


m' 


( e ) 


(5-44) 


On the other hand, the assumption that the wavelength of the ion acoustic 
wave and the electron thermal velocity are much greater than the Debye 
length and the phase velocity of the wave, respectively, i.e., 
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(5*45) 


( e )2 
% 


» k 2 s^ e ) 


» oi v 


.(e) 2 


gives the dispersion relation for an ion acoustic wave as 



( 2 ) 

The variation of ' with respect to 0^ gives the equation for 

the conservation of action, 


3t 3z 


= 0 


> 


N, 


~ I a / 
" 2 ? 


/ z* z\ 

Z„ 1 

3 / z* z\ 

( e k e 0 V*k ) C -°-> 

v, N, = - — 
k k 2 

dk \k e 0 K k e k) 
z ' ' 


j + c . c . , ( 5 . 47 ) 


2 

where the action, N , and the group velocity, v , are defined in 

K K 

Equation (3-13) • Explicit expressions for N, and vf can be obtained 

K JK 

by substituting Equation (5-42) into (5.47). 

5*5*2 The Coupled Mode Equations 
( 3) 

By using in Equation ( 5 . 15 ) and the first order equa- 

tion of motion [the first expression in Equation (5.4l)], we can write 
Equation (4.10) as 


iN, 


(It + \ fl) “k - (g '<$*’) 


e h 
k > 


(5.48) 


k,k,k 


where the first summation is over all permutations of k, k ', k" ; the 

( Qj) 

second is over the electron and ion species, and Q> „ is given by 

KK K 


(a) 2 

(a) n a H 
\kV ( ) e 0 


4s 


(a Vk" 

z z 


^ H/ 0 2 - 3k % 2s f, a) ) \ + 4i 2 -i s S a 7 


(5.49) 
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where (-) a is 1 for electrons and -1 for ions. Equation ( 5 . 48) can 
be written in the form of Equation (4.l4) with 



(5.50) 


c = Q (a) ) • 

zzz y. \C£ -psi/ 
-P, s ,i 


( 5 . 51 ) 


Equation ( 5 . 5 O) has been obtained from Equations (5.42) and (5-47). 

The role of the ordinary mode can also be considered in the non- 
linear interactions, as pointed out at the beginning of this section. 
Equation (4.l4) is therefore the coupled mode equation for any combina- 
tion of the electron plasma wave (p), the ion acoustic wave (i), and 
the ordinary wave (o) . The combinations which satisfy the synchronism 
conditions are (i) pump P-signal and idler P, I (see Figure 5-2), and 
(ii) pump O-signal and idler P, I (see Figure 5* 1 )- 

Case (i) ( C ppj) : Using Equation (5.43) for the electron plasma wave 

and Equation (5 .45) for the ion acoustic wave together with their linear 
dispersion relations given by Equations (5-44) and (5.46), the quantities 
I and C in Equations ( 5 . 50 ) and ( 5 . 51) can be written as 

K 


i: = 2ie, 


co 


k.O 

TO 2 “ 


N 



(5.52) 


C PPI _1 


e o q 


TeT z. 

3m ' k.O) 

1 N 


1 / 2 z, z\ 

f-T- r + 12 k k , 
(e)2 \ D P s/ ' 


where k is defined by 
D J 


-1 


k^ = (Debye length) 


(5-53) 


(5-54) 
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(b) 


Figure 5.2 Synchronism parallelogram for interaction of two plasma 

waves (pump and signal) and one ion acoustic wave (idler). 
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In deriving Equation ( 5 . 53 ), the synchronism condition of Equation (4.13) 
has been used. For the special case where the Debye length is much 
shorter than the wavelengths of the ion acoustic waves (pump and signal), 
the coupling coefficient in Equation (5-53) becomes that obtained earlier 

■ v 88 

by Tsytovich. 

Case (ii)(C ): In the limit of k Z - 0 for the ordinary wave, 

v v OPI' P 

the wave vectors of the signal and idler waves are equal and opposite 
[see Figure 5 . 1 (a)], i.e., 

k = k Z = k Z . (5-55) 

s 1 


By using the approximations and dispersion relations expressed by 
Equations (5.43) - (5.46) for the electron plasma and ion acoustic 
waves, the quantities I and C in Equations ( 5 . 5 O) and ( 5 . 51 ) can 
be reduced to 




(a) 2 


■S 


°‘S‘ 

k ,0 


(5.56) 


c -tfo! k: 
OPI 5 


D 1 

k 2. (aJJa ) 2 • 
a p,o 


(5-57) 


The coupling coefficient in Equation (5-57) has been obtained by other 

89 90 

methods by Goldman and Harker. 


5 .6 Application of the Lagrangian by Method II 

The purpose of the following linear theory is to demonstrate that 
the second order Lagrangian obtained by the second method [Equations 
(5*33) and (5-34)] is indeed the correct one. To do so, we shall prove 
that the linear dispersion relations obtained from the hydrodynamic 
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Lagrangian are the same, within the low temperature approximation as 
those obtained by the microscopic treatment of Chapter III. For 
simplicity, only the electron contribution will be considered in the 
linear calculation. A similar demonstration can be made for the non- 
linear theory. Since the procedure is the same as for linear 
theory, the demonstration for an electron plasma is omitted, and we 
shall consider instead the interaction of two electron plasma waves 
propagating nearly parallel to the static magnetic field with an Alfven 
wave propagating precisely parallel. For this interaction the ion 
motions will be included, since they play the dominant role in determining 
Alfven wave properties . 

5.6.1 Linear Theory 

As in Section 5 . 5.1, the Lagrangian can be written as 
Equation (5.40) for the linear theory, with A^ 2 ^ as in Equation (5-33)* 
The first order variables to be independently varied are 
and a , with the gauge cp = 0. After some albegraic manipulation, the 
Euler-Lagrange equation from the variation in yields 
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Cartesian coordinates (x,y,z) were used in Equation (5.35)* For sim_ 
plicity, the quantities k, a^ (hence e^ and b^) , and W in Equa- 

tion (5.58) are expressed in the rotating coordinates (R, L, Z), via 
the transformation of Equation (3.31). It should be noted that the 
transformation for 3/dW yields 


3/3W T = £ U" 1 a/9w. , (I = R,L,Z j = x,y,z) 

I , IJ J 


( 5 - 59 ) 


where U T 1 is the inverse of U_ . . 

Ij Id 


The variation of in Equation (5.4o) with respect to 


gives 


(03 03 + k s,+k s )? 

v k.O k,=Fl 11 z ||'* 
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k ,0 k,Tl' z || aw z m\k z k / ’ 


2 .,2 2 _ S JZ 


k . f \,_ k 


( oi +k s +k s ) £ = (cd + 03 )k s 5 — + (u) ^403 , )k s ^ — 

''HkjO 1 1 z || y ^k v k,0 k, 1 ' R i a«r R v k,0 k,-l' L 1 3w 


a? 

+ 203 „k s vrp + - \ . 
k,0 z || dW z m k 


(5.60) 


The Euler-Lagrange equation resulting from the a ^-variation yields 
the wave equation, as in Equation (3.19)> 


H X (n X e,) + k_ 


~k 


e^ = 0 , 
~k ’ 


.ij 


1 4 - 


ie -0) e 

0 k k 


(5.61) 


where the current density is given by 
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= - 1 'l( d . + s „ hq} = - k R( u z + s „ 1 ^) • 


(5-62) 


Equations (5-58), (5.60) and ( 5 . 61 ) are the required first order hydro- 
dynamic equations. 

If the mean value of the Larmor radius and the thermal velocity 
of the electrons are very small compared to the wavelength and the 
phase velocity of the waves, i.e., 


2 2 

kfs, « <0 
IX c 


k 2 s 
z II 


« 03 0) 

k,p k,v 


(p,v = 0,±1, or ±2) , ( 5 . 63 ) 


for a wave whose frequency is far from the cyclotron frequency 

» 03 c or « 03 c ), the second relation in Equation ( 5 . 63 ) ma Y 
be replaced by 



2 2 
k s « 03^ 
z II k ,0 


(5-64) 


Elimination of V,„ 5 from Equations ( 5 . 58 ) and ( 5 . 60 ) gives the equa- 

2 2 2 

tion of motion to first order in k /cn and k s /oi 03. , as 

x c z || nc,|x k,v 7 


^k ^k' * 


(5.65) 
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where the elements of the matrix are given by 
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With the use of Equations (5*58) and (5.64), the equivalent plasma 
permittivity matrix can be expressed as 


(5.66) 
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( 5 - 67 ) 


The wave dispersion relations are then obtained by substitution of Equa- 
tion (5.67) into Equation (5,61). 

The equivalent permittivity in Equation (5.67) can also be obtained 
from Equation ( 3 . 44 ), with the low temperature approximations of Equa- 
tion ( 5 . 63 ), by using the following series representation of the J 

n 

Bessel function: 
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( 5 . 68 ) 


J n U) 


A\ n v 

2 j mi r(nri-n+l) 

m=0 


Thus, when the low temperature approximation is made, the dispersion 
relations from the microscopic treatment reduce to those obtained in 
this section. 

5.6.2 Interaction of Electron Plasma and Alfven Waves 

The interaction mechanism of two plasma waves and an Alfven 
wave is as follows. In the magnetosonic regime, the Alfven wave pro- 
pagating parallel to the static magnetic field is a linearly polarized 
transverse wave.^ 1 If its electric field is polarized in the direction 
of the x-axis, this electric field can interact with the x-component 
of the electric field of an electron plasma wave propagating nearly 
parallel to the magnetic field to produce rf space charge variations in 
the x-direction. This rf space charge may in turn excite another electron 
plasma wave also propagating nearly parallel to the magnetic field. 

These electron plasma waves have the same azimuthal angle which may, 
of course, be chosen arbitrarily. Since the strongest interaction ob- 
viously occurs when the electric field components of these electron 
plasma waves normal to the magnetic field are aligned along the x-axis, 
we shall assume that the electron plasma waves propagate in the xz-plane 
as shown in Figure 5 > 3 . Therefore, all vectors are in either x, y or 
z-directions . For albegraic simplicity, drift velocities are ignored 
in the following. 

Linear Equations: First, for the Alfven wave, we have 


a) « a/ 1 ) 
k c 




( 5 . 69 ) 
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Figure 5.3 Synchronism parallelograms for interaction 
of two plasma waves (pump and signal) and 
one Alfven wave (idler). 
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To lowest order in 



transformation from (R,L,Z) to (x, y, z) 


coordinates reduces Equations (5-58) and ( 5 . 65 ), respectively, to 



(5.70) 



Similar equations for the ions are obtained from Equations (5.70) and 
( 5 . 71 ) by the procedure described in Section 5-4. Including the ion 


motion in Equations (5.61) and (5.67) gives the dispersion relation as 



(5-72) 
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Next, tor the electron plasma wave, we have the relations 
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To first order in k /k , Equation (5*58) becomes 
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(5.74) 


and Equation ( 5 . 65 ) becomes 
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Substitution of Equations (5.67) and ( 5 - 73 ) reduces the wave equation 
[Equation (5.61)] to 


xx x , xz z . 

V \ % * 0 ’ 




where 


are given by 
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In Equation (5«77)> the ion contributions have disappeared because 
[Equation (5.43)] . 

Coupled Mode Equations: In order to obtain the coupled mode equa- 

( 3) 

tions, the general expression for A £"/.// in Equation (5-34) has to be 

KK K 

reduced according to the conditions of the present case. Let the Alfv^n 
wave be the idler. Equation (5*69) for the Alfven waveband Equation 
(5-73) for the electron plasma wave then reduce , X^q^ and X^^ 

in Equation (5-35) to the forms 
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(5-78) 


Since substitution of Equations (5-70)> (5-7l)> (5*74) and (5*75) into 
Equation (5.78) gives a rather complicated expression, an additional 
simplifying assumption will be made that 


w 


» 00 


(i) 


(k,k' = electron plasma waves). 


(5.79) 


After eliminating e by use of Equation (5.76) , the coupled mode equa- 

K 

tion can be written in the form of Equation (4.l4) with 
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C PPA 


V 1 /s Q z ( a ) 

Z-f -P s 
■p,s \ 



(5.80) 


are given in Equation (5-77) with k = 


P> and Q kk^ 


are defined by 
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Here (-) a is 1 for electrons and -1 for ions. It should be noted 

that Equation ( 5 . 81 ) has two kinds of terms, one of order O(ccP) and 

_0 ( 0 ^ \ 

the other of order 0(<n ) . Since cju^ ' » the former terms for 

the ion equation can be neglected compared to those for the electron 

equation. Therefore, for £ the summation need be taken only 

a 

for the terms of order 0 (cd“^) . It should also be noted that the 

x-components of the electric fields of the electron plasma waves, and 

hence the coupling coefficient, C__„, vanish when k -* 0. 

’ PPA x 


5.7 Discussion 

In this chapter, two methods of formulation of the hydrodynamic 

Lagrangian have been presented and applied to several specific cases 

of wave-wave interactions. Lagrangians previously used by others for 

42,69,70 

the MHD approximation include only the second moment equation. 

Since our hydrodynamic Lagrangians include Maxwell's equations and the 
second moment equation (first method), or three moment equations 
(second method), they are more general. Fundamental equations for 
the MHD approximation can be obtained from our Lagrangian in Section 

Q o 

5.2.3 for spherical compression by introducing Ohm's Law. 

The second method in this chapter has very important advantages 
over the first. In particular, as demonstrated in Section 5*6, it 
does not involve the mathematical difficulties associated with the 


expansion of the Lagrangian in the first method, which effectively 
restrict us to the three special cases treated in Section 5*2. 

Another significant advantage is that in the case of nonlinear inter- 
action it can be applied when all interacting waves do not have the 
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same value of y. This situation cannot be handled by the first 
Lagrangian method. For example, although the small-signal proper- 
ties of electron plasma waves and compressional Alfv^n waves in an 
isotropic plasma can be studied by the first method using two differ- 
ent second order Lagrangians [Equations (5.15) and (5.21)], for 
interactions between them we must use only one Lagrangian; since 
these two kinds of waves have different Lagrangians for third 
order and above, a satisfactory higher order Lagrangian cannot be 
chosen. In the second Lagrangian method, however, this difficulty 
does not occur; the interaction of any combination of waves can be 
treated provided that the basic assumption that the wave phase velo- 
cities are greater than the electron thermal velocity is valid. 

Finally, we note that the second Lagrangian method includes all 
necessary moment equations, whereas the first contains only the second 
moment equation (the momentum equation), and must employ the other 
moment equations as subsidiary conditions. 
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Chapter VI 


CONCLUSIONS 


6.1 Review of the Research 

The main purpose of our work has been to develop a general theory 
of the averaged Lagrangian method to describe warm, collisionless 
plasmas, and to demonstrate the effectiveness of the method in analyzing 
wave-wave interactions. This has involved formulation of appropriate 
Lagrangians, and extension of previous work on the averaged Lagrangian 
method for cold plasmas. Treatment of specific cases has demonstrated 
the power of the method in analyzing wave-wave interactions by con- 
firming some results obtained previously by other methods, and pro- 
viding coupling coefficients for a number of new interactions. Since 
our developments cover both the microscopic treatment and the hydro- 
dynamic approximation to it, the averaged Lagrangian method is now 
available for use in all of the usual plasma descriptions: microscopic, 
hydrodynamic, and cold plasma. 

Among the advantages of the averaged Lagrangian method, as compared 
with the iterative method, are first, as demonstrated in Section 4.2, 
that the coupled mode equations can be obtained in canonical form 
from c^-variation of the Lagrangian without specifying the detailed 
form of the Lagrangian for* a particular system. In the iterative 
method, on the other hand, the coupled mode equations can be obtained 
only after detailed and lengthy algebra specific to the case under 
study. In the Lagrangian method, tedious algebra is required only in 
obtaining the coupling coefficients. Second, the Manley-Rowe rela- 
tions can be obtained by variation with respect to phase, and are 
automatically consistent with the coupled mode equations. 
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The second formulation of the hydrodynamic Lagrangian proposed 
in Section 5-3 has significant advatanges over the iterative method 
(here we mean direct use of the moment equations and Maxwell's equa- 
tions). As described at the beginning of Section 5 . 2 , expansion of 
the moment equations involves mathematical difficulties. Consequently, 
even for small-signal propagation, most of the work by direct use of 
the moment equations and Maxwell's equations has been limited to 

either the three cases in Section 5.2, or the special case of a 

92 

plasma in a strong static magnetic field, and the expansion has been 
carried out using an adiabatic equation of state of the form of Equation 
(C.l). For cases more general than those just mentioned, the validity 
of the adiabatic equation of state is in doubt. For example, an elec- 
trostatic wave has been shown [Section 5.2.l] to satisfy Equation 
(C.l) with Y = 3- As soon as its propagation direction departs from 
the z-direction, the assumptions in Equation ( 5 . 11) are violated, and 
Equation (5.12) is no longer justified. As demonstrated in Sections 
5*3 and 5-6, and discussed in Section 5.7, such difficulties are 
avoided in the second formulation of the hydrodynamic Lagrangian. 

The averaged method has been used to obtain the coupled mode equa- 
tions for all possible interactions of waves propagating nearly 
parallel, or precisely perpendicular to the static magnetic field, 
using the microscopic equations. Several interactions have also been 
treated in the hydrodynamic approximation. The new interactions 
that we have studied by the microscopic treatment are those among 
three circularly-polarized waves (except that among three right-hand 
circularly polarized waves); among two plasma waves and one circularly- 
polarized wave, and among one longitudinal cyclotron harmonic wave 
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and two ordinary cyclotron harmonic waves (pump || - signal and idler 
0,0) . In the hydrodynamic approximation, the only new interaction is 
that among two electron plasma waves and one Alfv^n wave. Since one 
of the principal aims of our applications of the method has been to 
compare the results, and the ease of obtaining them, with previous 
analyses by other methods, we have not carried out numerical calcula- 
tions. We will simply speculate here on some future directions that 
appear to be worth following. 

It is evident that rf space charge plays an important role in 
mode coupling. In order for nonlinear interaction to occur, at least 
one of the participating waves must have nonzero rf space charge. This 
being so, we may divide all interactions into four cases: those among 

three longitudinal waves; among two longitudinal and one transverse 
wave; among one longitudinal and two transverse waves, and among three 
quasiperpendicular waves which have nonzero axial components. Previous 
analyses have shown that relatively strong interactions occur among 
three longitudinal waves (e.g., plasma waves or longitudinal cyclotron 
harmonic waves). For interaction of three cyclotron harmonic waves, 

for instance, Harker and Crawford calculated the nonlinear growth rate 
due to a pump wave of 0.4 volt to be order of 1 - 10 db/cm over the 

range of cyclotron frequencies for a typical laboratory plasma of 

10 x xp 

n^ = 1.25 X 10 cm , when the electron temperature is 4 volts. The 

interactions of this category studied in our work are the P-P,P and 
II — 1 | , || in the microscopic treatment, and the P-P,I for the hydrodynamic 
approximation. In the P-P,I interaction we note that the range of 
frequency for synchronism is narrow, and the electron plasma waves 
and the ion acoustic wave have very low group velocity, so that plasma 
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inhomogeneity or collisional effects my invalidate our results. 

The next important interactions should be those involving two 
longitudinal waves and one transverse wave. Examples treated here 
are the R-P,P; L-P,P; P-R,P; and P-L,P for the microscopic description, 
and the O-P P and P-P,A in the hydrodynamic approximation. The mechanism 
of interaction by the microscopic treatment may be considered to be 
similar to that of interaction among three longitudinal waves, since 
longitudinal components of circularly-polarized waves can produce rf 
space charge in the direction of propagation when they propagate 
slightly oblique to the static magnetic field. In the case of pre- 
cisely parallel propagation the longitudinal components vanish, and 
so does the interaction, as shown in Table 4.1. Because these longi- 
tudinal components are small, we may expect these interactions to be 
weaker in general than those involving three longitudinal waves. How- 
ever, as mentioned in Section 4.5, the whistler and plasma waves can 
be unstable for a non-Maxwellian plasma, and can also carry positive 
or negative energy. It would consequently be of great interest to make 
further detailed analyses of the R-P,P and P-R,P cases when the R-wave 
represents a whistler. The character of the 0-P, P interaction is 
essentially the same as that of the P-P,P interaction since, in the 
dipole approximation, the ordinary wave behaves like a plasma wave 
of infinite wavelength. The interaction may be expected to be corres- 
pondingly strong. In fact, this interaction is currently considered 
to be an important process in ionospheric modification experiments 

Q 

with high power sounders , and in producing plasmas by use of lasers. 

As shown in Figure 5*3( a )> the P-P,A interaction occurs due to coupling 
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of the components of electric field perpendicular to the static mag- 
netic field. The interaction consequently vanishes for parallel 
propagation. Since the analysis in Section 5 -6 was for the Alfv^n 
wave within the magnetosonic regime (o^ « the range of fre- 

quencies of the electron plasma waves for synchronism is relatively 
narrow. In practice, slight inhomogeneity of the magnetic field and 
plasma would cause strong deviation from the theoretical predictions. 

Interactions considered in our work involving one longitudinal 
wave and two transverse waves are the R-R,P; L-L,P and P-P,L, for parallel 
propagation, and the 0— 1| , 0 and ||-0,0 for perpendicular propagation. One 
of the features of these interactions is that the second order electron 
gyrations due to two transverse waves ^ X are cou P led to the 

rf space charge of the longitudinal wave. Therefore, contrary to the 
case of interactions involving two longitudinal waves, interactions 
can occur among waves propagating precisely parallel to the magnetic 
field. Since only one longitudinal wave is involved, the interaction 
may be expected to be very weak. 

The final case of interaction is among circularly-polarized waves 
propagating slightly oblique to the static magnetic field, whose 
electric field components in the direction of wave propagation produce 
small amounts of rf space charge. Examples are the R-R,R; R-R,L and 
L-R,L. The rf space charge is much weaker than for the first two, and 
hence no strong interactions can be expected. In the ionosphere, the 
growth rates due to the R-R,P; L-L,P or P-R,L are of order 10 ~^~ 10 ~ ^ db/m 
for a pump wave field strength of e p ~ 0.1 m volt/m, ^ and that due to 

interaction among three whistlers is of order 10 - ^ db/m or less for a 

27 

pump field of e^ ~ 1 m volt/m. Previous calculations have, however, 


been performed assuming cold plasma, or by use of moment theory. It 
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should be remembered that electron plasma waves and whistlers can be 
unstable, and can also carry either positive or negative energy when 
the electron velocity distribution is nonMaxwellian. Further detailed 
studies are consequently necessary for interactions involving these 
waves under magnetospheric conditions. 

6.2 Suggestions for Future Work 

The next important step following our work is to obtain numerical 
solutions for the growth rates for comparison with observations made 
in the laboratory and the ionosphere. Most previous numerical 
solutions have been obtained for cold plasmas, or within the hydro- 
dynamic approximation. Those obtained by the microscopic treatment 
for warm plasmas are potentially of greater interest because the micro- 
scopic description generally offers a broader range of possibilities 
for synchronism. Also, interactions involving longitudinal waves or 
whistlers, which would otherwise be weak, may be much more important 
in plasmas having nonMaxwellian particle velocity distributions. 

A large number of interactions which can be treated in a straight- 
forward manner by the averaged Lagrangian method still remain to be 
studied. In the microscopic treatment, for example, the interactions 
between waves propagating parallel to the static magnetic field, and 
those propagating perpendicular should be examined. In the hydro- 
dynamic approximation, there are numerous combinations of electromag- 
netic waves and ion waves which are likely to couple strongly, and to 
be observable experimentally. 

In our work, the averaged Lagrangian method has been developed for 
coherent three-wave interaction in an infinite, collisionless warm 
plasma in a uniform static magnetic field.. Some extensions of the 
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present theory to less restricted conditions would be of interest: 


1. Higher Order Wave-Wave Interactions : Waves which do not 

satisfy the synchronism conditions for three-wave interaction can still 

be coupled by higher order wave-wave interactions in which more than 

three waves interact. Even when three-wave interaction occurs among 

existing modes, such higher order wave-wave interactions may not be 

negligible, since they may have growth rates comparable to those of three- 

wave interaction. To describe n-wave interaction, the Lagrangian must 

th 

be expanded up to n order in perturbation. 

2. Background Reaction : If the reaction on the assumed homogeneous 

and time- invariant plasma is taken into account, the Euler-Lagrange 

equations for the perturbed state should be obtainable from variations 

of the expanded Lagrangian with respect to r(= r + ?), 9(=$ + cp ) and 

— * — -0 *<- 1 

43 

A(= Aq + a ) , instead of J* , cp^ and a of Section 2.2. The Euler- 

Lagrange equations from Z 2 should then yield equations equivalent to 

4 

those of quasilinear theory. An alternative way to take account of the 
reaction on the background is to assume the equilibrium velocity distri- 
bution function, f Q , to vary slowly in time and space, and to take 

variations in and a . This method is being developed at 

96 

Stanford by Galloway. 

3. Nonlinear Wave-Particle Interaction ; Next in order of complexity 
is nonlinear wave-particle interaction. This corresponds to the beating 
of two waves to form a sum or difference wave which then interacts with 
particles in resonance with it. As for wave-wave interaction, the 
Lagrangian must be expanded beyond second order in perturbation to take 
such effects into account. The main difference from wave-wave inter- 
action is that the particle velocity distribution itself is modified by 
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the waves. Consequently, the technique developed to include the back- 
ground reaction has to be extended . 

4. Nonun iform Plasmas : The expanded Lagrangians in Section 2.2 

are applicable to inhomogeneous plasmas. Homogeneity has not been 
assumed in formulating the hydrodynamic Lagrangian by the second method 
in Section 5.3. Therefore, microscopic or macroscopic properties of 
inhomogeneous plasmas can be described by using these Lagrangians. 

Among the subjects of interest is stability analysis. Suydam obtained a 

97 

stability criterion for pinch in the MHD approximation. Similar 
criteria for microscopic and macroscopic plasmas could be derived by 
obtaining the Hamiltonians from the Lagrangians, and by finding the 
conditions for minimizing the potential energies. It is noted, in passing, 
that, since the field quantities are functions of r^ + £ , the 

98 

Lagrangian in Equation (2.15) is in the form of a correlation function. 
Interesting properties of the Fourier transform of the correlation 
function would suggest the possibility of simpler analysis than the 
direct use of Maxwell's equations and the Vlasov equation for the 
microscopic treatment or the moment equations for the hydrodynamic 
approximation . 

5. Bounded Plasmas : Amongst the complications not treated compre- 

99 

hensively is that of plasma boundaries. So far, Laval et al, and Spithas 
and Manheimer 100 , have considered the three-wave interaction of three 
coherent plasma waves in a waveguide with an infinite magnetic field, 
while P6rulli 101 has made an extensive study of such processes from 
both a theoretical and experimental viewpoint. All of these studies 
have used the conventional iterative method. The averaged Lagrangian 
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method has been employed in recent Stanford work by Larsen for the 
interaction of three coherent plasma waves in an infinitely long cold 
plasma column with and without an infinite static magnetic field. The 
three-wave interaction work should be extended to the case of a plasma 
column of finite length in a finite magnetic field. 

6. Random Phase Approximation : So far, we have discussed only 

interaction of coherent waves. As mentioned in Section 1 . 2 , however, 
the random phase approximation must be used for turbulent plasmas. For 
randomly phased waves, only quantities averaged over a statistical 
ensemble are meaningful. Since simple averaging makes the LHS of 
Equation ( 4 . 10 ) vanish, the coupled mode equations cannot be used. 
Instead, the wave kinetic equation, which describes the rate of change 
of the plasmon distribution function or quantum density, is used. 1 ^ 

The wave kinetic equation is the ensemble average of the action transfer 
equation [Equation (4.9)3*^^ and the plasmon distribution function or 
quantum density is the action, N . Therefore, the coupling coefficients 

XV 

105 

obtained in this work lead directly to those for turbulent plasmas. 
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Appendix A 


DERIVATION OF EQUATION 2.18 


Figure A.l shows two pairs of corresponding cells defined in 
Section 2.2. One is at time t, and the other is the same pair after 
the infinitesimal time At. The equilibrium cell and nonequilibrium 
cell are located at (r , v ) and at (r,v), respectively, at time t, 
and at (r',v') and (r / ,v / ) at time t ' - t + At . The displacement 
vectors are ( ?, £) at t, and (F 7 , F 7 ) t 7 . 

The displacement vectors at t and t 7 are given by 



(A.l) 


The position vectors r / and r 7 at t 7 can be written as 

= £ + Olo + ^ = lo + ^ At > ( 

where use has been made of v = v^ + E. Substituting Equation (A. 2) 
into Equation (A.l), and subtracting the first expression in Equation 
(A.l) from the second, gives 


A? = ?' -•§ = 5 At . (A. 3) 

rvi ro ~ <-N -/ 

The velocity displacement vector, E, is therefore the time rate of change 
of the displacement vector The time rate of change of §, on the 

other hand, can be found by the Taylor expansion of § 7 about §, and 
substitution of ^ from § 7 . It yields 

% - [S + • 7 > k + (t ■ ’v)i] • (A- 1 *) 

Equation (2.18) results from Equations (A. 3) and (A. 4) . 
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Figure A.l Time evolution of a corresponding cell pair. 
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Appendix B 


INTEGRAL REPRESENTATION OF THE COUPLING COEFFICIENTS 
FOR PERPENDICULARLY PROPAGATING WAVES 


In this appendix we shall derive expressions for the coupling co- 
efficients in Equations (4. 32) -( 4 . 34) in convenient forms for computa- 
tion. This will be done by expressing the Q' s in Equation (4.28) in 
terms of integrals. To obtain integral expressions for the DG' s and G' s 
we must solve the equation of motion expressed by Equation (3.18) or 


(3.32). 

In the absence of 



solving Equation (3-32) for ^ gives 



°k,0 °k 




(B.l) 


where 


„RR 


* LL* q exp{-i[o^ i(f<j) k )-k i v i s i n( ( j)- ( t)k)] / cb c } 

= D , (J , = — — 1 : — 7 — 7 * 


r» p _ 

k,0 k “ ^-kjO^-k “ 2im co c sin( jta^/co c ) 


^k +ir _ r - k x v j_l [" — _ 1 

X J d<t)|exp(-i<j>) + exp|^i(cD k <j) + k^sin^/coj , 


„ _,RL * LR* q 

\,o G k = D -k,o G -k = “ TH5 


k x V ±exP { i ( <M> k ) "i [<\( <M> k ) -k x v x sin( <j)-(J) k ) ] /oj c } 

a3 c CO k Sin ( ,T( \ / “c) 


- < l )_( t ) k + Tr 


r r T k T]l _ 

I d(j) exp[ i(cc^()H-k v sin <►)/ CD ] 

A 1 1 


(B.2) 




Equation (B.2) in turn gives 
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RR LL* q 
G k = G _ k = TJi 


exp {-i[cu k ( (j)-(t) k ) -k^sin( (J)-({) k ) ] /cd c } 

4 si An<\/%) 


^k +rt _ _ +3T =r = k v 1 = = 

x J d<j) exp(i<|>)J_ d<() Jexp(-i<|>) - exp[ ito^-k^sin^/o^] , 






„RL LR* 

G, = G 
k -k 


q k 1 v _ L exp ^ 2 i< i > k -i t^( _ $ k ) - k J _ v J _sin( ]/co.} 


Bin 


sin 2 (rt^ k /a3 c ) 


/ 

J *\-« 


M* _ _ 

exp(i(j)) I d(j) exp[ i ( k v sin(j))/<D c ] . 
■?-n 


(B.3) 


R L Z II J. Z 

The unitary transformation from (e k , e k , e k ) to (e k , e k , e k ) ^ described in 


Section 3 . 5 .I, requires that 


a " - I (of + of + G f ex »(- 21 i) + °f e ^( 21 M • 

G k" * - \ (of - of + exp( -2i6 k ) - G“‘exp(2i^ k )) . (B.4) 


Equation (B.4) also applies to D G . Substitution of Equations 

K j U K 

(B.2) and (B. 3 ) into Equation (B.4) gives 
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D G »" = 
kj Ok 2im 


exp (-i (cD k ( <j>-<}> k ) -k^sin(({i-(l) k ) ] /<»<.} 

tJD c sin ( ir % /ai c ) 


(b-m 

/ 

*V<L 


"V 11 _ 


cos(<|>— <j)-<j> k ) exp[i(cn k (|> + k^sin^/a^] , 


... exp{-i[a^((l>-<}) k )-k i V x sin((j)-<|> k )]/cD c } 

n __ . — — 

k, 0 k - 2im (jD c sin(rt<xi k /oo c ) 


* 

4m 


T Tk +n _ _ _ 

I d(j) sin( <J>— <J>— 4> k ) exp[ 1 ( 0^4 + k^sin^/coj , 

■'K-* 

exp £-i [ 0 D k ( 4 )-<l> k )-k i v x sin (<t>-<i> k )]/cD c } 


2 2 

a) Sin ( % ) 


_ Jh-* = _ _ _ _ 

I d<|) I d<t> cos((j)-4) exp [ i ( co^fj) - k^v^sin^AuJ , 

Jt i Jl 


J k 4m 


4-4 k “« 

q exp {- i [ ( 4-4 k ) -k x v a s in ( 4-4 k ) ] /a3 c } 


2 2 

a) c sin (jta^/cD.) 


_ rk* = 

d<j) I d(j) sin( exp[ 1(0^4 - k v sin(j>)/a3 c ] . 


Equation (4.28) is then written as 

2 


<kV' - - 35 8 o k l T ( T kV *M-i[<V*- k ? i si ”(H I '') / “c’> > 

CO 


c 

2 


Q kk'k" = B S oK (t(2o) c /co k ,)S k - (l-V'/o^')^]) > 

CO 


where S and T are defined by 
K K 


(B.5) 


(B.6) 
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S. = 


ex p(- i “ k < l ) k /“ c ) ^ k+Tr 
sin( I ta^/co c r J 

< P- ( P k -« 


d(j) cos (<M>-<i> k ") exp{-i[a^<j> + k^v^sin<J)]/a> c } 


exp(-iai <j), /oo ) ^k +Tt _ _ _ _ 

k = l 2 >* * - £ _ d( j> d< t)cos(H 

sin («V“c) >* 


' k -<t> k /' ) exp [-■ 1 [ £ \^“ k j. v i sin ^] /uy c 5' 

(B.7) 


Substitution of Equations (B .6) and (B.7) into Equations 
an integral representation for the coupling coefficients 
and , as desired. 

oo„ 


(4. 32 ) -(4. 34) gives 
C || II II ’ C „oo^ C 0,|0^ 
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Appendix C 


EXPANSION OF THE STATE EQUATIONS 


The purpose of this appendix is to obtain the elements of the 
tensor *s second and third order in terms of the perturbation J|, 

/■n_> 

for each of the three illustrative cases in Section 5.2. For all 
three, the state equation [Equations (5- 1 3)> (5-l8) and (5.23)] is of 
the form 


d 

dt 



= 0 , 


r = y-i , 


(C.l) 


where y is the adiabatic compression constant. In order to expand 
Equation (C.l), we should first transform the total time derivative 
d/dt, in the (r,t) coordinates to that in the (^^q) coordinates. 
Equation ( 5 . 8 ) immediately gives 


a • a : a_ ^ <L_ 

at Sfi'ar ' Sr 


Therefore, it follows that 


(C.2) 


d 

dt 



(C.3) 


Equation (C.3) states that the rate of change of a function along the 
perturbed trajectory is the same as that along the equilibrium tra- 


jectory. 

From Equation ( 5 . 8 ) it follows that 



1 ? a § 
t - 1 + 7 -i + 2 [<’•!)■ - 7 i ;7 y + ^ 


' lay x )’ 


(cA) 


where the subscript 0 on x,y,z has been dropped. In order to expand 
Equation (C.l), let 
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s = s 0 + s l + s 2 + " ' * $ = 1 + ^1 + ^2 + ’ ■ ' ‘ ( c -5) 

Here C can be expressed in terms of E from Equation (C. 4 ). Substi- 
n ***» 

tuting Equations (C. 4 ) and (C.5) into Equation (C.l), and expanding it, 
gives the s^ as 

s i = " s o c i » s 2 = ' s o (C 2 _c i ) > s 3 = - s 0 [C 3 -Ci( 2 C 2 -Ci)],.-. • (c. 6 ) 


(i) Parallel Compression : From Equation (5.I3), T = 0 and 2 for 

s and s , respectively. By Equation (5 . 11 ), the Jacobian in Equation 
(C. 4 ) reduces to 


i = 1 + 


5 z 


(c.7) 


Equation (C.6) then becomes 


S jJL “ S x 2 " S x 3 


= 0 , 


s ,2 = 3s »o (sr) ’ s ,| 3 = - hs ,o(4rf > (c - 8) 


s . = -2s 
||! || 0 Sz 


(ii) Cylindrical Compression : From Equation (5.18), T = 1 and 0 for 

and s^ , respectively. By the assumptions in Equation (5.16), the 
Jacobian for this case becomes 


$ = 1 + 7 


55 55 55 55 

P + _JL _ _EE _Z 

~ 5x 5y 5y 5x 


(C.9) 


Equation (C.6) then yields 
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(c .10) 


P as d? x 1 

S xl = _S lO (7 ± * S) J s x 2 = S xoL (7 x ’ ~ " ax“ r dy~ + 5y" aFJ 


s = -s 
±3 J-0 


r « 3? as as as ~\ 

<Vi> [<Vi> 2 - 2 ^ i/ +2 5F 


s =s =s = ...=0 
111 II 2 II 3 


(iii) Spherical Compression : From Equation (5-23 ) > r = 2/3- 

Equations (C.4) - (C.6), the s' s are obtained as 


s„ = - — s„ 
1 3 0 


7 -l> s 2 “ 1 s o[f (7 ' I) 2 + 7 l ;7 i] > 


- 3 s o[lr < v -s> 3 - § < 7 • - 2 1 • (1 * I)] 


From 


(C .11) 
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Nonlinear Interaction Between Circularly Polarized Waves and Electron Plasma Waves* 

H. Kim, K. J. Harker, and F. W. Crawford 
Institute for Plasma Research, Stanford University, Stanford, California 94305 
(Received 6 August 1970) 


This paper presents a general analysis of nonlinear three-wave interaction in a hot magnetoplasma 
between circularly polarized waves propagating at small angles to the static magnetic field, and electron 
plasma (Langmuir) waves propagating parallel. First, the coupled mode equations are derived by iterative 
solution of the Vlasov equation. Simplified expressions for the coupling coefficients are then obtained 
for all of the wave combinations for which the frequency and wave number synchronism conditions can 
be satisfied, and the wave coupling coefficients are nonzero. These comprise interactions among three 
right-hand polarized waves, two circularly polarized waves and one Langmuir wave, and three Langmuir 
waves. The paper concludes with a brief discussion of the specific cases most likely to be worth subjecting 
to detailed numerical analysis. 


I. INTRODUCTION 

The nonlinear interaction of three or more waves in a 
magnetoplasma has been a subject of extensive in- 
vestigations in recent years . 1 Initially, for simplicity, 
most of these investigations employed either cold 
plasma theory, or macroscopic theory based on moments 
of the Vlasov equation. More recently, interest has 
progressed to studies in which the Vlasov equation itself 
is used. The primary reason for this is the fact that even 
certain important small-signal phenomena, such as 
collisionless Landau and cyclotron damping, and the 
propagation of cyclotron harmonic waves, only emerge 
from a microscopic treatment. It is consequently 
essential to begin studies of nonlinear interactions 
involving these phenomena from the Vlasov equation. 
An additional reason is that certain new effects appear 
which have no counterpart in cold-plasma theory. In 
particular, for nonMaxwellian charged particle velocity 
distributions, the possibility of nonlinear interaction 
between waves carrying positive and negative small- 
signal energies exists. This can lead to explosive 
instability , 2 3 in which all of the wave components grow 
at the expense of the charged particle energy. 


Nonlinear plasma wave propagation is a notoriously 
complicated subject because of the wide variety of 
phenomena that may occur 4 : Individual waves may 
suffer nonlinear damping due to wave-particle interac- 
tion, and two or more waves may interact with the 
particles, or with each other in pure wave-wave interac 
tion. In each case, for the interaction to be strong, 
certain synchronism conditions must be satisfied 
between such quantities as the wave frequencies, wave 
numbers, cyclotron frequency, and charged particle 
velocities. It seems likely that progress in such a com- 
plex situation will be made most easily by studying 
each effect separately in its region of greatest strength, 
and comparing the results for a given set of plasma 
parameters to determine which interaction effects 
should dominate. The final predictions, and if feasible 
the separate analyses leading to them, should then be 
subjected to experimental verification. 

In this scheme, the authors have opted to concentrate 
on wave-wave interactions, in which the charged 
particle velocity distribution is assumed to remain 
steady in time. To further restrict the field, attention is 
limited to the high-frequency propagation branches 
obtained ignoring ion motions, and to propagation at 
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6 = 1 r/2 or 8 x 0 , i.e., corresponding to interaction 
between principal modes propagating either perpendic- 
ular or parallel to the static magnetic field. So far, the 
case of perpendicular propagation has been examined 
in most detail by the authors and others. Theoretical 
studies have been published of three-wave interactions 
between cold plasma, transverse (extraordinary) 
waves, 5 6 and between hot plasma, longitudinal (cyclo- 
tron harmonic) waves, 7 ' 8 and have received encouraging 
support from laboratory experiments. 910 For approxi- 
mately parallel propagation, a cold-plasma treatment 
has been given for the interaction of right-hand polarized 
transverse waves, 11 and the interaction between 
circularly polarized modes and a longitudinal electron 
plasma (Langmuir) wave has been studied, using 
moment theory for the longitudinal wave. 1213 The 
purpose of the present paper is to extend this work to 
include hot-plasma effects stemming from the Vlasov 
equation. 

The general procedure for analyzing three-wave 
interactions adopted here is as follows. Maxwell’s 
equations and the Vlasov equation are first expanded to 
first order. This yields equations for the field quantities 
and a small-signal dispersion relation. The expansion is 
then taken to second order, and all second-order terms 
are approximated by substituting first-order quantities 
in them. This can be done for three waves provided 
that the small-signal dispersion relations allow their 
frequencies and wave numbers to satisfy certain 
synchronism conditions. This procedure is carried 


through in Sec. II for a general charged particle velocity 
distribution, and leads to a set of coupled mode equa- 
tions. 

In Sec. Ill, these coupled mode equations are 
specialized to three cases. The first is that of three right- 
hand polarized waves propagating almost parallel to the 
static magnetic field (there is no interaction in the limit 
of parallel propagation). Three left-hand polarized 
waves need not be considered since the topology of the 
dispersion characteristics does not allow the synchronism 
conditions to be satisfied. Mixtures of left- and right- 
hand polarization do not couple. The second is that of 
two circularly polarized waves interacting with a 
Langmuir wave, all three waves propagating strictly 
parallel to the magnetic field. The synchronism con- 
ditions can be satisfied for three different combinations 
of wave types, and these are treated individually. The 
case of one circularly polarized wave and two Langmuir 
waves propagating parallel is of no interest as coupling 
cannot occur. The third case considered is that of 
interaction among three Langmuir waves, and our 
treatment confirms the form of the relevant coupling 
coefficient given recently by Dysthe. 14 

Having obtained the general coupling coefficients 
and coupled mode equations, specific problems of 
mixing and parametric amplification may be attacked. 
Although no detailed numerical solutions are given in 
this paper, some comments on the procedure to obtain 
them are given in Sec. T1I.D, and the likely results are 
discussed in Sec. IV. 


II. THEORY 
A. Basic Equations 

Our system is described by Maxwell’s equations 

vxe=-aH/a/, vxh= g+(ds/dT), (l) 

supplemented by an expression for the convection current density 

3 = — jFVrfV. (2) 

The evolution of the charged particle velocity distribution function F is described by the Vlasov equation 

( dF/dt ) + V- VF- ( £+ V x H) • ( dF/dV ) -U c {dF/d<t>) =0, (3) 

where we have used cylindrical coordinates in velocity space defined such that 

V X =V cos<t>, V p =-Vsm<t>, (4) 

and the static magnetic field Bo is oriented along the z axis. 

Equations (l)-(4) have been written in terms of the normalized variables 

£=—eE/mcu 0 , H=— eB/mwo, d= — J/noec, ,V = v/c, 

X=wqx/c, T=u 0 t, S2 = or/ wo, U c =~eBo/moio, K=ke/o>o, (5) 

where w 0 is the electron plasma frequency, and the remaining symbols have their usual meanings. 

Taking the spatial-temporal Fourier transforms of Eqs. (1) and (3) yields 

K v x £ T =S2 T H T , K t xH 7 =— 3j-jn y £ y , 

7n T F 7 -y(v-K 7 )F T +fi c (aF T /a<#.)-(£ T -t-vxH T )-(aF„/av)= £ £ a +VxH a ).(«ydV), (6) 

a,0 
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where it has been assumed that the synchronism conditions 

a y =Q a +U 0 , K 7 =K a +Kjs, (7) 

for three-wave interaction are satisfied. The first two expressions of Eq. (6) are linear. The third has linear terms on 
the LHS, and nonlinear driving terms on the RHS. 

In what follows, we shall restrict our attention to circularly polarized waves propagating along or at small angles 
to the static magnetic field, and electron plasma waves propagating along it. Under these conditions, we may write 
the wave number components as {KJ) xy i x , K y 9 vy i„, K y i 2 ), where i,, i„, i„ are unit vectors. If we introduce the com- 
plex transverse vectors 

S T = S IT +ys„, H y = H iy +jH vy , 3y = 3ty+jdvy= ~ jF y V 2 dVdV z exp(-i<#>), (81 

appropriate to propagation of right-hand polarized waves, then to lowest order in 9 Eq. (6) yields the expressions 

(K y ya y *) (£ iy 9 y -& y ) +£ y = ( j/Oy)g y , {K y y 2iy) (w*+&*e y ) +& iy =j$jo, y , (9) 

and , 

dF 

(Q y -K y V,)F y +jQ c —? 


exp( -j<t>) +S 7 exp#] [ F, 




+j |s„ [<J+ ^ [9 y exp j<b+9* exp(-#)] [v ^ - V 

= ~ j S { S2 “ [ a+ k C0 “ cxpjW exp( ^ )] ( Fs kr~ V 5fJ] 

+i(exp/'<f>) (z c ,b a +je,J a ^j +i[exp(— #) ] (t^*b a -jZ-*d a E s +K VK y ) [0 V exp #+0 Y * exp( -#)]F r . 

(10) 

Here, F ( , is the unperturbed electron velocity distribution function, and 0., has been written for 9 zy +j0 uy . The 
operators a and b„, and the velocity-dependent quantity d a , are defined by 

a=d/dV s , b a = (d/dV) + ( KJ9 a ) [V(d/dV,) - VAd/dV ) ], 4= F-'[1-(FX,/U.)]- (11) 

B. Iterative Solution 

Equation ( 10) may be solved iteratively to any desired order of approximation. The first iterate is obtained by 
setting Fp and F y to zero on the RHS and solving to obtain 

(12) 

(13) 

This result can be used as an improved guess for F s and F, on the RHS of Eq. (10). After performing a total of 
three iterations, we obtain all terms quadratic in £ and linear in 9. We may then write F y as a power series in 8 as 

F T = (U r n+G T 2i) expj0+(G,i O +C T 2o) + (G : v i_ 1 +G 72 _i) exp (—#) + ( other terms), (14) 

where G y u, G yl0 , G t2 i, and G y20 are given by 


where T y „ is defined bv 


T yn = U y — K-, V z —nQ c . 


G y n j 


. [ b y e 


2T y i + 29. 


Kydy£> z y ^ / 

2 Q y T y0 T yl L 7 ° V 


V Z ~~V T 7 ~) +n y Va 


dV dVj 
Gy io = — ids*? {o./T y ^)-\-9 y *Z y {K y V/^T y xT y Q)b y ~\ F 0 , 


]} c, 


, = -!(£«£-/) ( 


jrhjr (h+dft) — b a + y9y (b'.+d a )T 0 _ l %+ {b 0 +2d 0 ) V b a 

1 yl 1 70 al J- 7O J- 7I l ali a2 


_ V’ 

+ 7V ba T^T» 


2 KpOp 


.)f„- 


Ke«s,») ( T y r'b a T 0 0 ->oH- Tyr'aT'.r'b „) Fo, 


9gTyiT al 

^ 720 = t( 8 -^*£<i)C^ 70 _1 (^+ < W T al ~ l ba-\-T y (T l (ba-\-d a ) T 0 -{~ l bf\Fa—ZiJ2,iiT y ir'<i(T 0 <r x -\-T a < ) ~ v )aF< i . 


(IS) 
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Substituting Eq. (14) into Eqs. (8) and (9), we obtain 

(K y ya y *) (8 y s 27 -s 7 ) +e 7 = (-2 tj/Oy) ismvdVAGw+GTu) , 


{K 2 /2Q, 2 ) (£,«/+£-A)+£ !r = i-2rj/n y ) J jVdVV.dV,(G yW +G yW ) . (16) 

Substitution of the expressions defined in Eq. (15) finally reduces these to 

[1- (X T 2 AV) + (27rA2 T ) /F 2 rfW 1 (6E 0 /2r Tl )]S T +*.A&«Y= (.-2rj/O y ){{V i dVdV,G ya , 

W-fiy*^+^~2^RF 0 VdVdVJT^ n = (-lwj/%) nVdVVJVfir*, (17) 

where 'F T is given b) ; 

%= (K y 2 /Q y 2 ) -2rffF 0 VdVdVA {K y V s /T yl Q y 2 ) + (K y 2 vy2Q y )[_(T y0 2 T yl )~'+ (K y V z /a y T y0 T yl *)l). (18) 

C. Velocity-Space Integration 

Equation (18) is the desired general form of the coupled mode equations, describing how the electric field com- 
ponents and of the y wave are influenced by those of the a and 0 components. To solve a specific problem, it 
is necessary to specify the unperturbed velocity distribution F 0 and to carry out the integrations. Considerable 
progress can be made if the not unduly restrictive condition is imposed that F 0 is of the separable form 

Fo(V, V z )=F 0J .(V)F 0z (V s ), 2irjFo±(V) VdV= 1. (19) 

This will be assumed in what follows. 

It is also necessary to obtain expressions analogous to those of Eq. (17) for the a and 0 waves. The method of 
doing so is as follows. If we wish y to represent one of the lower frequency waves, rather than the highest as in 
Eq. (7) , it is only necessary to change the sign of 0 in the foregoing analysis. To demonstrate the procedure, and in 
anticipation of the requirements of Sec. Ill, we shall make this change. After carrying out appropriate velocity- 
space integrations, Eq. (17) then reduces to 

D/Z !y +’i> y 0 y *& y = Q Sy e, a 0p*+R y e, Z aS ! p*. (20) 

The quantities D y and D y r , and the coupling coefficients P y , Q y p, and R y , are given by 
Z) T = (AV/fV) - (1 /a y yfdV 3 F 0 l(T y0 /T yl ) +K y 2 ( P/27V)], 

Z5/=2[l-/</E 2 E 03 (l/r 70 2 )], 


/\ = 


8Q a QpQ y 


X) Kpdp J 


[ 4T T d 

dVJPo. \ ~ -£ - £ ■ +K a KpK y V 4 — L(T al T gi T yl )-'(T a r'+T e ,-«)] 


+ 


2K y 2 V 2 (T. 


T a \T$\Ty\ \ 


+ p°) +2*7 P 

1 al * y\J 


Ty0 (T a r'+T s r')+ (TfH-'+Ttr 1 ) 




TpoTg\T y \ 


+ 


2KpK y Y 2 

T a \Tg\T y \ 


Tao T y i ^ 2 7'„o 2Tgi 

\_Tg 0 Tp o T a i T, 


J 


'7&~ 


/ f , r/ P 2 / K y T a o FaT-,t)\ 

2a,0, J ‘ J |_7V V T ai T yl ) 


-K a K y V 2 — (T al Tp 0 T y ,)-' 
oV , 


R y =-2jfdV ! F< H (d/dV z ) ( T a0 Tp 0 T y0 ,)-'. (21) 

The coupling coefficient Qp y is obtained from Q y p by interchanging 0 and y. In the above expressions, a bar denotes 
an average value of the quantity concerned, i.e. , of V 2 and V 4 . As may readily be guessed, Eqs. (20) and (21) re- 
sult from some very lengthy manipulations. Attention is drawn to some of the more important points in the deriva- 
tion in an Appendix to the paper. 

The significance of D y and D y p is worth noting. If we consider small-signal waves propagating parallel to the 
static magnetic field, the expressions in Eq. (20) reduce to Z> 7 6 7 ~0FyZ?/'E 72 , and we recognize the coefficients as 
the well-known hot-plasma dispersion relations for circularly polarized transverse waves, and longitudinal electron 
plasma (Langmuir) waves, respectively. 
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III. APPLICATIONS 

We are now ready to apply the general results of 
Eq. (20) to specific cases. Three possibilities im- 
mediately come to mind: first, that all three waves 
might be circularly polarized transverse waves; second, 
that two circularly polarized waves might interact with 
a Langmuir wave; and third, that three Langmuir 
waves might interact with each other. We will indicate 
the appropriate coupling coefficients for these interac- 
tions separately, and then discuss how they may be 
used in the solution of mixing or parametric amplifica- 
tion problems in inhomogeneous plasmas. To emphasize 
the type of wave concerned, we will introduce the 
variables 

£>a R = £ > xa-\~j&ua, — £ 2 a ’j&yay 

£/=S 2 „, («>0), (22) 

referring to right- and left-hand circularly polarized 
waves and Langmuir waves, respectively. 

A. Three Circularly Polarized Waves 

Eliminating £, from Eq. (20), and neglecting terms 
higher than first order in 6 and second order in £ yields 

£ 7 %* = C„ K £„% fi *, (23) 

where the coupling coefficient Crr is given by 
C RR =P y - £ (WAs P ) 

y.0 

¥0=1— sifr'SFtodv, 

x[r„-H-(jyt*/2 r*r*)(r*-'+7y->)]. (24) 

It is interesting to note that Crr —* 0 as 0„, 0 y —> 0. Thus, 

there is no interaction for purely parallel propagation, 


Table I. Acceptable combinations of waves for 
interaction (0=0). 



Signal and 

Pump (p) 

idler ( s , i ) 

R 

R,P 

L 

L,P 

P 

R,L 


and this property is independent of the precise form of 
the velocity distribution. 

In previous work, 11 the authors have considered 
interaction in a cold plasma. It is possible to retrieve the 
expression for C RR derived there by evaluating Eq. (24) 
for the delta-function velocity distribution E 0 (F, V z ) = 
8(V)5(V Z ) /2ttV appropriate to a cold plasma. 

B. Two Circularly Polarized Waves and a 
Langmuir Wave 

In this case, interaction is possible for 9 = 0. Equation 
(20) then simplifies to 

D y Z y =Q.JlJ-,/*, £>/£/= QjtSoV. (25) 

Of the several possible combinations of left- and right- 
hand polarized waves ( L , R) and Langmuir waves (P), 
there are only the three given in Table I for which the 
synchronism conditions can be satisfied, and the coupling 
coefficients are nonzero. 13 For illustration purposes, we 
have called the highest frequency component the 
“pump” wave, and the lower frequencies the “signal” 
and “idler,” analogous to conventional parametric 
amplification terminology. 15 We will now consider the 
acceptable combinations separately. 


1. Pump R — Signal and Idler R, P (or P, R) 

The required coupled mode equations are obtained from Eq. (25) as 

D/Z* = - C RP Z «£, p , D,'% k = CwS^fi,*, DfZf = C RP Z **£„", ( 26) 


where D, R , D P R , and DP are given by D y in Eq. (21) after replacing y by s, p, and i, respectively. The coupling 
coefficient C R p is given bv 


J 

2Q.SI, 


f nr P r 2 / K,T p0 K p T t0 \ - 3 

J 7 W P * dV z 


(Tpi T l0 P,,i 


(27) 


The first and second expressions in Eq. (26) have been obtained from the first expression of Eq. (25) by replacing 
(a, 0, y) with (s, — i, p) and ( p , i, s), respectively. The last expression in Eq. (26) derives from the second of 
• Eq. (25) by replacing (a, 0, y) with (p, s, i). 


2. Pump L— Signal and Idler L, P (or P, L) 

The results for this case can be obtained simply by changing the sign of K v , K„ K u fi p , I2„ and 12, in Eq. (26). 
When this is done, the electric field vectors change from £„ B , £/, £, p , to £/*, £, L *, £, p *, respectively. Upon making 
these substitutions, and then taking the complex conjugate, Eq. (26) becomes 


D p L e p L =-C LP £, t L &i p , 


DpZp=C L ,ZP*?> v L , 


D i p & i p =c LP e„ L *s p L , 


(28) 
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where D-, L is defined as 

ZV'= 1 - (K y 2 /Qy 2 ) - (1/a;) SdVJ'ol(T y0 /T y ^)+K y H (y=p, s), (29) 

and the coupling coefficient Clp is given by 

C L p=(j/ 2 n p Q s )fdV l F 0 A( 2 /T i o i )l(K s T p0 /T p - l )-(K p T, 0 /T l - l )l-K p K,V*(d/dV s )(T p - i T, 0 T,- l )-'\. (30) 


3. Pump P — Signal and Idler R, L (or L, R) 


The results for this case are obtained from Eq. (26) by changing the sign of ft, and K,, and then making the 
interchange (p, s, i, p). When this is done, the field vectors change from £,*, S P K , 8 < p , to S, L *, 8„ R , S p p , 

respectively. Upon making these substitutions, we obtain from Eq. (26) the coupled mode equations 

D/£/=-C RL e,PZ s R , £>/£.« = A R e, L =C flz .S**S/, (31) 


where the coupling coefficient Cr L is given by 

2 


Crl — 


3 (, IVF [ 2 (K.T, KiT. 0 \ 
2Q&J ‘ ;F<ij Lt 5 o 2 \ T,-i + T tl ) 


+K s K i V i — (T^T^T^)-' 

ov z 


(32) 


C. Three Langmuir Waves 

For this case, we put 00 = 0 = 0, in Eq. (20) to obtain 

D/G/=2?,G« P S0 P *. (33) 

The coupled mode equations take the form 

D/E,/= -Cpp& p G, p , D P S P =C PP £/S, P *, 

A P G, P =C PP G/G/*, (34) 

where D p p , D. p , and DP, are given by Eq. (21) with y 
replaced by p, s, and i, respectively, and the coupling 
coefficient Cpp is given by 

Cpp= -IjfdVfoPd/dV,) (T p0 T, 0 T i(l )-'. (35) 

Here, the expressions in Eq. (34) have been obtained by 
replacing (a, 0, y) in Eq. (33) with ($, —i,p), (p, i, s), 
and (p, s, i), respectively. 

D. Procedure for Solution of Problems 

The coupling coefficients derived in Eqs. (24), (27), 
(30), (32), and (35) constitute the primary informa- 
tion needed for detailed numerical studies of the non- 
linear interactions involved. The picture is only' com- 
pleted, however, when one inserts in them wave 
numbers and frequencies satisfying the synchronism 
conditions and the small-signal dispersion relations. 
These have been given elsewhere for the cold-plasma 
approximation to the results of Sec. III. A, 11 and the 
hydrodynamic approximation to those of Sec. III. B. 13 
They would have to be recalculated for any other 
velocity distribution chosen. 

If we assume that the synchronism problem has been 
solved at some point in a weakly inhomogeneous 
plasma for which the dispersion relation D y is of the 
form 

D y =D-,lK(X),n, X], (36) 

then the most general problem is to determine how 8, 


varies. It will be noted that all of our wave equations 
are in the general form 

D y & y = C,G„G0. (37) 

Using the technique of Berk and Book, 16 we can expand 
these equations around the point of noninteraction to 




Fig. 1 . Synchronism conditions for nonlinear interaction 
between three right-hand polarized waves in a cold plasma. 
(Heavy portions of curves indicate regions over which the pump 
and idler can vary. Arrows indicate signal and idler variation as 
pump increases.) (a) Pump in upper branch; (b) Pump in lower 
branch. 
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Fig. 2. Synchronism conditions for nonlinear interaction between two circularly polarized waves (R or L ) , and a plasma wave ( P ) 
in a cold plasma. QCut-off frequencies flt/.z,— ( l-|-ficV4) 1/J ±tie/2.] (a) Pump R — Signal and idler R, R; (b) Pump /. -Signal and 
idler L, P; (c) Pump P — Signal and idler R, L. 


obtain 

d&ydDy k r d_ (dDy\ &Dy aKr"] dZy dDy 

dXdK y + * e ' y Lax WJ + dKS dX J dT ao v 

= -/C Y S„£ s . (38) 

Thus, once the dispersion relations and coupling co- 
efficients are known, one may determine through Eq. 
(38) the variation of each wave amplitude in time and 
space, in homogeneous or weakly inhomogeneous 
plasmas, and for whatever type of problem is of interest, 
e.g., wave-wave mixing, or parametric amplification. 
Several examples of the procedure are given in Ref. 13. 

IV. DISCUSSION 

A very important point to note concerning the 
coupled mode equations [Eqs. (23), (26), (28), (31), 
and (34)] derived for the various cases discussed in 
Sec. Ill is that the coupling coefficients for the signal 
and idler equations are invariably equal in magnitude 
and sign, and equal in magnitude but opposite in sign to 
that for the pump equation. This is sufficient to 
guarantee that the Manley-Rowe relations are sat- 
isfied. 15 The applicability of these relations to the 
problems studied here is not only of theoretical im- 
portance, but also constitutes a check on the accuracy 
of the extremely lengthy algebra and manipulations 
leading to the coupling coefficients. 

Although these coupling coefficients [Eqs. (24), (27), 
(30), (32), and (35)] appear to be quite complicated, 
they are in fact straightforward to program for com- 
puter solution. Numerical results can then be obtained 
for wave growth due to mixing, or parametric amplifica- 
tion, as soon as the detailed form of the unperturbed 
charged particle velocity distribution F 0 is specified. 
We have not yet carried out such a program, and will 


simply speculate here on some directions that appear to 
be worth following. We will do so by considering what 
additional features will appear when previous treat- 
ments are extended to include hot-plasma effects. 

The first consideration in three-wave interaction 
problems is to determine what opportunities the small- 
signal dispersion relations offer for the frequency and 
wave number synchronism conditions to be satisfied. 
We shall examine first the case of three right-hand 
polarized waves propagating in a cold plasma. There are 
then two situations in which synchronism can be 
achieved. 11 These are shown in Fig. 1. In the first, the 
pump and signal lie on the upper branch, and the idler 
lies on the lower (whistler) branch. In the second, all 
three waves lie on the lower branch. 

Suppose now that the plasma electron velocity dis- 
tribution is Maxwellian. Thermal effects will make only 
a very slight difference to the shape of the upper 
branch, and will introduce no collisionless damping 
there, since the phase velocity exceeds the velocity of 
light. The lower branch will be susceptible to collision- 
less damping, however, and this will become particularly 
strong for tt c - Sl<JfF ( . Here, V c («1 for our non- 
relativistic theory to be valid) is the thermal velocity 
normalized to the velocity of light. We may reasonably 
conclude that as V e increases from zero, the weak 
parametric growth predicted by the cold-plasma theory 
will be progressively quenched. The quenching effect 
will be greatest if any of the interaction frequencies 
approach so the case with all three waves on the 
lower branch will probably suffer most. 

The situation becomes much more interesting if the 
plasma electron velocity distribution is nonMaxwellian. 
The upper branch will still be almost unchanged from 
its cold-plasma shape, but the lower branch can be 
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strongly modified and become unstable. 17 The possibility 
exists of both positive and negative energy modes 
occurring, and these might lead to explosive instability. 
It would be of great interest to determine the relative 
growth rates of the small-signal and nonlinear in- 
stabilities to see which would dominate in experimentally 
realizable situations. 

So far, we have only discussed the interaction of three 
right-hand polarized waves. We now consider the three 
cases for which coupling coefficients were derived in 
Sec. III. B. Typical synchronism conditions for a cold 
plasma are shown in Fig. 2. As shown in Figs. 2(a) and 
(b), the combinations of R-R, P and L-L, P do not 
look very interesting for further study: even the 
Langmuir wave has a phase velocity of the order of the 
velocity of light, so the effects of nonrelativistic thermal 
motions will probably be small for all three wave 
components. 

The R-R, P case is actually richer in synchronism 
possibilities than has been indicated in Fig. 2(a). For 
the plasma frequency lower than the cyclotron fre- 
quency, i.e. , Cl c > 1, it is possible for the pump and signal 
to lie on the lower branch. For an isotropic Maxwellian 
velocity distribution, all three waves would then be 
subject to collisionless damping. For anisotropic dis- 
tributions, a wide variety of interesting combinations 
involving damped and growing waves would be possible, 
and would require extensive investigation. We may also 
remark that, although for cold plasma it is not possible 
to achieve synchronism with the pump and signal in the 


upper and lower branches, respectively, this condition 
could be achieved with suitable velocity distributions, 
and is open to investigation. 

The remaining combination in Fig. 2 is that of 
P-R, L. We note that the plasma wave has a phase 
velocity of the order of the velocity of light, so that it 
will be very little affected by the inclusion of thermal 
velocities. The whistler branch may be damped and/or 
unstable, depending on the velocity distribution chosen, 
but it seems unlikely that significant nonlinear effects 
will occur. It should, incidentally always be remembered 
that in cases such as this one, where one of the waves 
( P ) has an extremely low group velocity, plasma in- 
homogeneity will severely limit the distance over which 
the synchronism conditions required for interaction can 
be satisfied. 

The final case treated in Sec. Ill, i.e., that of three 
Langmuir Staves, has no solution in cold-plasma theory, 
since the synchronism conditions obviously cannot be 
satisfied in this limiting case of dispersionless plasma 
oscillations. For nonMaxwellian velocity distributions, 
a wide variety of damped and growing solutions can 
occur, however, and present a very important topic for 
further detailed nonlinear studies. To take a specific 
example, it would be valuable to examine the case of an 
approximately monoenergetic beam interacting with a 
Maxwellian plasma to determine to what extent the 
fastest-growing small-signal component is damped non- 
linearly by wave-wave interaction, as opposed to wave- 
particle interaction. 


APPENDIX: DERIVATION OF THE COUPLING COEFFICIENTS P y , Q y$ , Q 0y 

After performing double integrations with respect to V and V z , and changing the sign of 0, Eq. (17) reduces to 
the forms expressed by Eq. (20), but with coupling coefficients defined by 




^ ff dVdV z F a (KJ y b a T a c'(M-de) 


V 3 

TylTyl 


— KyBybpTpC 1 (ba—da) 


V 3 

Tyl Ty 0 


V V 2 

+K a 6 a b„ —— „ ■ (bg— 2 dg) — ■ 

a\*- a2 * yl 


V V 2 2Kg9 0 V 2 

7>i7> 0 K T yl + b ° T al T yX 


)■ 


Qyl>= ( irj/n y )f!dVdV z Flb a T a r'a( V 2 /T yX ) -aT„-'b.{ V 2 /T yl )\ 

Q» y = (*j/ 0 y) ndVdV z F 0 lb a T al -'(b e -d s ) ( V z V/T y 0 ) -b g T^-'{b a ~d a ) {V,V/T yl )]. 


(Al) 


By use of Eqs. (11) and (19), the expression for Q yS may be reduced to that of Eq. (21). It may also be shown 
that Qg y is related to Q y p simply by interchange of /3 and y. The main problem is that of reducing P y . To accomplish 
this, we note that the wave number synchronism condition implies that 


KJ> a — K-.0-,. 


(A2) 


It follows that the equation for P y can be separated into two parts: one involving K/flg, and the other involving 
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K-,6- y. After differentiation of the integrands, we have 

[QTaO 


jfr (?^T> ccr..r«)-- ( r„r„:- 3 ) 

+2K a W\T^T y ,r'+2K a Kf i VHT, 1 T m T, i )-'+2K a K.S^^-(T a r l - -2 KfV 2 ^°° (7V+7V 1 ) 

1 yl~ \ It 91/ 1 01* 00* yl 


~~ (T'.r'-Ter 1 ) 

1 70-* 7l 

Tp 0 


H-a^a,!* — {(r. 1 r 7 l )- , C(7' a 2 7’, 1 )- I -(7'm7', 1 )- , -(7' 41 r 7 „)-']!+2^y 5 v 

dF* 1 C,\ l -,\l yO 


+2K a K y V- 


Tgg 

TylTyO 




T yl ~'+ Ty 0 ~') + 2K a K^( T. A Ty 0 )-' 


T B r'-T a r')-2Kt, 2 V 2 


TaO 

2 til l T y\T yft 


+ 2A fl A 7 P 


l S\TylTyl 


( T.,r'+Ty ( r') +4 a 3 a y r- 


r al 7\ 


[(r 0 ,r,„)- | -(r m r Yl )- 1 - (z>,7\„)- 1 ] . (A3) 


P T is to be rearranged so that it shows symmetry in /3 and y. Each of the two integrands in Eq. (A3) has terms 
of zero, second, and fourth order in V. The symmetry in 7 and /3 of the zero and fourth-order terms is comparatively 
easy to show in the form of the first two terms in the large curly bracket of Eq. (21). The symmetric property of 
the second-order terms, i.e., those in V 2 , is not an easy matter to demonstrate, however. To do so, we note first that 
symmetry of the last terms in the integrands of Eq. (A3) is evident from comparison with the V* terms. By use 
of the synchronism conditions to eliminate K a , the rest of the 1- terms may be divided into three groups, involving 
As 2 , Ky 2 , and KaK y . After some tedious algebra, all of the l ' 2 terms then reduce to the symmetric form expressed 


in Eq. (21). 
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